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ABSTRACT 

Aims. We aim to provide atmospheric parameters and rotational velocities of a large sample of O- and early B-type 
stars, analysed in a homogeneous and consistent manner, for use in constraining theoretical models. 
Methods. Atmospheric parameters, stellar masses and rotational velocities have been estimated for approximately 250 
early B-type stars in the Large (LMC) and Small (SMC) Magellanic Clouds from high-resolution VLT-FLAMES data 
using the non-LTE TLUSTY model atmosphere code. This data set has been supplemented with our previous analyses 
of some 50 O-type stars (Mokiem et al. 2006, 2007) and 100 narrow-lined early B-type stars (Hunter et al. 2006, Trundle 
et al. 2007) from the same survey, providing a sample of ~400 early-type objects. 

Results. Comparison of the rotational velocities with evolutionary tracks suggest that the end of core hydrogen burning 
occurs later than currently predicted. We also show that the large number of the luminous blue supergiants observed in 
the fields are unlikely to have directly evolved from main-sequence massive O-type stars as neither their low rotational 
velocities or position on the H-R diagram are predicted. We suggest that blue-loops or mass-transfer binary systems may 
populate the blue supergiant regime. By comparing the rotational velocity distributions of the Magellanic Cloud stars to 
a similar Galactic sample we find that (at 3a confidence level) massive stars (above SMq) in the SMC rotate faster than 
those in the solar neighbourhood. However there appears to be no significant difference between the rotational velocity 
distributions in the Galaxy and the LMC. We find that the v sin i distributions in the SMC and LMC can modelled with 
an intrinsic rotational velocity distribution which is a Gaussian peaking at 175kms~^ (SMC) and lOOkms"^ (LMC) 
with a i half width of 150 km s^^. We find that in NGC346 in the SMC, the 10-25 Mq main-sequence stars appear to 
rotate faster than their higher mass counterparts. It is not expected that O-type stars spin down significantly through 
angular momentum loss via stellar winds at SMC metallicity, and hence this could be a refiection of mass dependent 
birth spin rates. Recently Yoon et al. (2006) have determined rates of GRBs by modelling rapidly rotating massive star 
progenitors. Our measured rotational velocity distribution for the 10-25Mq stars is peaked at slightly higher velocities 
than they assume, supporting the idea that GRBs could come from rapid rotators with initial masses as low as 14 Mq 
at low metallicities. 
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1. Introduction 

The evolution of massive stars has traditionally been de- 
scribed in terms of mass and metallicity. However in the 
past 20 years it has become evident that rotational effects 
are equally important, from star formation through to their 
deaths in supernovae explosions and associated gamma- 
ray bursts (Woosley & Heger f^DUS)) . For example, Heger 
& Langer (|2000p discuss how rotation changes the lifetime 
and evolution of massive stars on the Hertzsprung-Russell 
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diagram through rotationally induced mixing between the 
core and envelope. 

Rotational effects are also believed to be metallicity de- 
pendent. At low metallicity, line-driven winds are weaker 
and mass-loss rates lower (Kudritzki et al. 1987J Kudritzki 
& Puis [Mini Vink et al. lMITl Mokiem et al. 20GEI. Hence a 
star will lose less angular momentum and thereby maintain 
higher rotational velocities on the main-sequence. This pro- 
motes rotational mixing over a longer period of an object's 
evolution which in turn extends its lifetime. Additionally 
nucleosynthetically processed material will be mixed 
from the core into the photosphere and hence abundance 
anomalies are expected, with helium and nitrogen enhance- 
ments and carbon and oxygen depletions. Recent observa- 
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tions of massive stars have revealed abundance anomalies 
in 0-stars (Bouret et al. IMBl Mokiem et al. [Ml^ [Wf|) 
B-type supergiants (Dufton et al. |2005, Hunter et al. 120071 
Lennon et al. [20031 and Trundle fc Lenn on [2005]) . B-type 
main-sequence stars (Hunter et al. 120071 Korn et al. 120021 
I2005P , A-type supergiants (Venn lT999)) and lower mass stars 
(Smiljanic et al. [20n5|) . 

Maeder & Meynet (|200ip have included rotation and 
metallicity effects in stellar evolutionary models to explain 
both the relative frequency of blue and red supergiants as 
well as observed abundance trends in A-type supergiants at 
low metallicity. However much remains uncertain. For ex- 
ample, they adopt initial rotational velocities of 300 km s~^ , 
which equates to an average main-sequence rotational ve- 
locity of 220-240 km s^^ (there being little change in the 
rotational velocity over the main-sequence lifetime of a B- 
type star), but statistically significant observational sam- 
ples are necessary both to estimate typical rotational ve- 
locities as well as to calibrate the effect of metallicity on 
the rotational velocities. Additionally rotation has been de- 
scribed as the key ingredient in enabling single stars to 
explode with associated long-duration gamma-ray bursts 
by Woosley & Heger (|2006p . although rotational veloci- 
ties close to 400kms~^ are required and hence gamma- 
ray bursts would be relatively rare events compared to su- 
pernovae. The requirement that little angular momentum 
is lost during their stellar lifetime obviously favours low 
metallicity regimes with their smaller mass-loss rates. 

A large survey of massive stars has been undertaken at 
the European Southern Observatory using the Fibre Large 
Array Multi-Element Spectrograph (FLAMES) on the 8.2m 
Keuyen Very Large Telesope. The observations are dis- 
cussed in detail in Evans et al. (|2005l hereafter Paper I) 
and Evans et al. (|2006[ hereafter Paper II) and are sum- 
marised in Sect. [2l Briefly, some 750 stars have been ob- 
served with mainly O and early B spectral types in the 
Galaxy and the Large (LMC) and Small (SMC) Magellanic 
Clouds. The analysis of many of these objects is well under- 
way or completed. Mokiem et al. (|2006l I2007P have derived 
atmospheric parameters and helium abundances for the ma- 
jority of the 0-type stars in the sample and discuss the de- 
pendence of mass-loss rates on metallicity and calibrate the 
wind-momentum luminosity relation. Hunter et al. (|2007p 
and Trundle et al. (|2007p have derived atmospheric parame- 
ters and metal abundances (CNO in particular) for a subset 
of the narrow-lined B-type main-sequence and supergiant 
stars finding that significant mixing has occurred in many 
of the objects. Dufton et al. (|2006i hereafter Paper HI) have 
estimated atmospheric parameters and rotational velocities 
of the entire Galactic sample and report that more massive 
stars have lower rotational velocities (as expected from en- 
hanced mass-loss) and that cluster objects appear to rotate 
significantly faster than objects in the Galactic field. 

In order to examine the effects of rotation, here we 
perform a similar analysis to Paper HI and estimate at- 
mospheric parameters and rotational velocities for the 
Magellanic Cloud sample of B-type stars. The principle 
intention of performing such an analysis was to pro- 
vide rotational velocities for a large sample of early- 
type stars in a homogenous and consistent manner. 
This dataset can then be used to answer several 
questions of importance to massive star evolution: 



Where does the end of core hydrogen burning 
occur and how can large populations of post core 
hydrogen burning objects be explained? 
Are the rotational velocities of early B-type stars 
dependent on the metallicity of the host environ- 
ment? 

Are mass-loss effects important at low metallic- 
ity and how do these affect the predictions of 
gamma-ray burst rates from massive star pro- 
genitors? 

In Sect. [2] we briefly describe the FLAMES survey and 
our selection criteria for the objects discussed here. In 
Sect. 13] we present our methods for deriving the atmo- 
spheric parameters and projected rotational velocities. In 
Sect. [4] the rotational velocity distributions of the 
different subsets of the sample (for example, bina- 
ries, Be- type stars, high mass stars and supergiants) 
are compared. In Sect. [5] the relevant subsets are 
utilised to answer the three questions posed above. 
Finally in Sect. [6] we present the principal conclusions from 
our analysis. 

2. Observations 

As part of a European Southern Observatory (ESO) Large 
Programme on the Very Large Telescope (VLT, Paranal, 
Chile) observations were obtained for over 700 stars us- 
ing the Fibre Large Array Multi-Element Spectrograph 
(FLAMES), with the majority of these objects being O- or 
early B-type stars. Observations were taken towards three 
Galactic clusters, NGC 6611, NGC 3293 and NGC 4755; two 
LMC clusters. Nil and NGC 2004 and two SMC clusters, 
NGC 346 and NGC 330. Details of the observational fields, 
target selection and data reduction were discussed in Paper 
I and Paper iQ- The FLAMES spectroscopy covered the 
wavelength regions 3850-4755A and 6380-6620A and was 
supplemented with additional observations from the Fiber- 
fed Extended Range Optical Spectrograph (FEROS) and 
the Ultraviolet and Visual Echelle Spectrograph (UVES), 
both of which had a more extensive wavelength coverage. 
To maintain consistency throughout the analysis we only 
consider the FEROS and UVES data in the same wave- 
length range as that of the FLAMES observations. 

2.1. Radial velocity corrections and binarity 

Observations were taken in six wavelength settings in order 
to achieve our desired wavelength coverage at high spec- 
tral resolution, with multiple exposures (typically six) be- 
ing taken in each setting. However, in cases where ex- 
posures are of significantly lower quality than the 
norm, these have been excluded. 

To increase the signal to noise (S/N) ratio of the data, it 
is possible to combine the individual exposures in each set- 

^ We note that in Paper II the radial distance of the NGC 346 
objects from the centre of the association was calculated based 
on the NGC 346 centre obtained from the SIMBAD database 
(operated at the CDS, Strasbourg, France) which is actually the 
centre of the ionised shell of gas and is offset from the centre of 
the OB association. In Table[T](only available online at the CDS) 
the radial distances are listed with the centre taken as the posi- 
tion of object 435 from the Massey et al. (I1989|l catalogue of OB 
stars in NGC 346; a(2000) = 00 59 04.49 5(2000) = -72 10 24.7. 
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ting. However, it is first necessary to check for radial veloc- 
ity shifts between exposures, an indicator of binarity. IDL 
procedures have been developed to cross-correlate each ex- 
posure in the FLAMES data and identify significant radial 
velocity shifts at the 3a level (see Hunter et al. I2007p . We 
note that this procedure is superior to that used in Papers 

I and H where binaries were identified by eye from the com- 
bined spectra of the exposures in each wavelength setting. 
Our cross-correlation procedure has identified additional bi- 
naries not previously detected in Paper II. Additionally for 
several stars we find that the radial velocity shifts detected 
in Paper II are not significant at the 3cr level and hence we 
do not consider them as binary stars. We note that most of 
these objects were identified as possible binaries in Paper 

II or had Be classifications and large rotational velocities, 
which makes binary detection uncertain. After correction 
for velocity shifts, IRAF routines were used to combine 
the exposures and simultaneously remove cosmic rays. This 
procedure was not possible for the double lined spectro- 
scopic binaries and the analysis of these systems (where 
this was possible) is discussed in Sect. 13.5] Throughout this 
paper all radial velocity variables arc considered as binary 
objects although our temporal resolution is not sufficient 
to derive periods and other orbital parameters to confirm 
if they truly are all binary stars. 

2.2. Selection criteria 

We have selected objects with spectral types later than 
O 9 and earlier than B 9 which can be analysed with 
our non-LTE TLUSTY model atmosphere grid which has 
an effective temperature range from 12 000 K to 35 000 K 
(see Sect. 13.11) . The majority of the 0-type stellar sam- 
ple observed in the FLAMES survey has been analysed by 
Mokiem et al. p006l [2007|) for the SMC and LMC 0-type 
stars respectively. 

We have attempted to derive parameters for every ob- 
ject in our spectral type range but this was not possible 
for a small minority. Typically objects were excluded if it 
was not possible to separate individual spectra in a dou- 
ble lined binary system, if the spectra were not of sufficient 
quality to reliably estimate velocity corrections in the case 
of single lined spectroscopic binaries or if the stellar spectra 
appeared to be variable between exposures. Additionally a 
few objects have been excluded as shell absorption contam- 
inated both the hydrogen and helium lines and no metal 
lines were well enough observed to allow for a reliable de- 
termination of the rotational velocity. In Table [2] we list 
those objects within our spectral type range that have been 
excluded. These objects form less than 6% of the sample, 
with the majority being binaries and having either Be clas- 
sifications or uncertain spectral types. 



2.3. Field populations 

In Paper II images of the four regions are displayed and it 
is clear that the majority of our objects do not lie within 
the core of the clusters, due to constraints in positioning 
the fibres of the FLAMES instrument. Specifically for the 
NGC 330 and NGC 2004 fields, our samples are probably 
dominated by field objects. NGC 346 does not have a com- 
pact core and it is difficult to determine the extent of the 
cluster. Indeed, NGC 346 is often described as an OB as- 



Table 2. Objects with spectral types later than 9 and 
earlier than B 9 observed in the FLAMES survey that have 
been excluded from this analysis. 



Cluster 



Excluded stars 



NGC 346 017, 038, 060, 086, 105, 111, 115 

NGC 330 030, 077, 088, 092, 093, 115, 117 

Nil 005,030,053,067,099,112 

NGC 2004 019, 028, 033, 034, 037, 038, 072 



sociation rather than as a bound cluster (see for example 
Gouliermis et al. I2006p . N 11 is dominated by two distinct 
regions, LH9 and LHIO, where sequential star formation 
has occurred. Additionally there are many smaller pock- 
ets of star formation within this region complicating the 
picture further. 

The supergiant objects in NGC 330 are not located 
close to or in the core of NGC 330; they instead appear 
to be randomly distributed in the field and may indi- 
cate that they are indeed field stars. However, it should 
be noted that almost no supergiants are observed towards 
NGC 346, implying an age difference between NGC 346 and 
NGC 330. Similarly to NGC 330 the supergiants observed 
towards Nil and NGC 2004 are not clustered in one re- 
gion. Additionally in Nil the 0-type stars are not clus- 
tered in the younger LH 10 region, indicating that the OB 
association is probably not bound. 

As such our four Magellanic Cloud samples are probably 
better described as unbound associations with the associa- 
tion 'evaporating' as their OB stars become part of the field 
population. Sirianni et al. (|2002p have obtained colours for 
the stars in both the central and outer regions of NGC 330 
and find remarkable similarity between the regions imply- 
ing that NGC 330 may extend well beyond its dense core. 
The recent review of Bastian & Gieles (|2006p describes the 
various phases of cluster disruption, in particular that of 
'infant mortality', where through gas expulsion during clus- 
ter formation the cluster becomes unbound, and if the star 
formation efficiency is low, the cluster can become unbound 
on the order of a few tens of Myr (see the recent work of 
Goodwin & Bastian 120061 and references therein) . 



3. Analysis 

3.1. Non-LTE tlusty model atmosphere grids 

The non-LTE tlusty model atmosphere grids (see Hubeny 
& Lanz 119951 and references therein) used in this analysis 
have previously been discussed - see Hunter et al. (|2007p 
for an overview of the grids and Dufton et al. (|2005p for a 
detailed discussiorfl. 

Briefly, model atmosphere grids have been calculated 
for the analysis of B-type stellar spectra, covering the ef- 
fective temperature range from 12 000 K to 35 000 K in steps 
of not more than 1500K, a surface gravity range from 
the Eddington limit up to 4.5 dex, in steps of 0.15 dex 
and microturbulence values of 0, 5, 10, 20 and 30kms^^. 
Four model atmosphere grids have been generated with 
metallicities representative of the Galaxy ([Fe/H]=7.5 dex). 



See also http://star.pst. qub.ac.uk| 
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the LMC (7.2 dex), the SMC (6.8 dex) and a lower metal- 
licity regime (6.4 dex). A nominal helium abundance of 
11.0 dex has been adopted throughout. Theoretical spectra 
and metal line equivalent widths have then been calculated 
for each grid point. 

3.2. Atmospheric parameters 

Hydrogen and helium lines are observable in the spectra of 
all our targets and these lines have been used as the pri- 
mary indicators of effective temperature and surface grav- 
ity. IDL procedures have been developed which fit model 
spectra (convolved with an appropriate stellar broadening 
function) of the appropriate metallicity to the observed hy- 
drogen and helium lines and calculate the region of best 
fit (by chi-squared minimisation) in effective temperature - 
surface gravity space. As these parameters are interdepen- 
dent it is necessary to use an iterative method. In Paper 
III equivalent width measurements of the hydrogen and he- 
lium lines were used to constrain the atmospheric param- 
eters but we believe that the 'profile fitting' methodology 
we have adopted here is an improvement as it uses both 
information on the line strength and shape. We estimate 
that this methodology results in a typical uncertainty of 
1 500 K for the effective temperature, with that for the sur- 
face gravity being discussed below. 

The He I spectrum is a suitable effective temperature 
indicator at spectral types of approximately B 3 and later, 
while the He ii spectrum can be used to constrain the tem- 
perature of stars with a spectral type of B 1 or earlier. Hence 
there remains a range of spectral types where the He ii 
lines are not observed and the He I spectra are not sensitive 
to temperature. For spectra with narrow metal absorption 
lines, where two ionization stages of silicon were observed, 
it is possible to use the silicon ionization balance to con- 
strain the effective temperature, with the surface gravity 
again being estimated from the hydrogen lines. Hunter et 
al. (12007 1 and Trundle et al. (|2007p have presented analyses 
of the narrow lined B-type stars in the FLAMES survey and 
their atmospheric parameters have been adopted. For stars 
with large projected rotational velocities in the B 1-3 spec- 
tral type range, it was not possible to observe two stages of 
silicon ionization. As in Paper III, we have adopted effective 
temperature-spectral type calibrations. These are based on 
the analysis of the narrow lined objects and are presented 
in Trundle et al. for both LMC and SMC metallicities. 

Both the H(5 and H7 absorption lines have been ob- 
served in our spectra and these have been used to estimate 
the surface gravity by the method of profile fitting described 
above. In the majority of cases we find excellent agreement 
between the estimate derived from each line with differences 
typically being less that 0.15 dex. The effective temperature 
and surface gravity estimates are directly correlated and an 
over estimate of the effective temperature by 1 000 K will 
lead to an over estimate of the surface gravity by approx- 
imately 0.1 dex. Hence we believe that our surface gravity 
estimates should be typically accurate to 0.2-0.25 dex. 

It has previously been found that there is a correlation 
between surface gravity and microturbulence with super- 
giants typically having greater values of microturbulence 
than giant and main-sequence stars (see for example Gies 
& Lambert 1992). Representative microturbulence values of 
10 and 5kms~^ have been adopted for supergiants and the 
less evolved objects respectively. It should be noted that 



this choice has a negligible effect on the atmospheric pa- 
rameters and only affects the projected rotational velocity 
estimate when the measured values are less than 20kms~"'^. 
The atmospheric parameters are listed in Tables [3]|6l only 
available online at the CDS. In these tables we give the 
star identifier and spectral type (from Paper 11), the ef- 
fective temperature, surface gravity, projected rotational 
velocity, adopted methodology, luminosity, mass and addi- 
tionally indicate if the star is a radial velocity variable and 
hence probably in a binary system. 

3.3. Projected rotational velocities 

Several methods are available to estimate projected rota- 
tional velocities (wsini). The most commonly used in the 
case of OB-type stars are based on a direct measure of the 
FWHM of the spectral Hu es (viz . Slettebak et aL [T9751 Abt 
et al. ,2002; Herr ero et al. HMIl Strom et al. I^UUS)) or pro- 
file fitting (Gra yfT976l Paper HI), while Penny (fT996l) and 
Howarth et al. (jl997) applied a cross correlation method 
to lUE spectra. Recently, Simon-Diaz & Herrero ( 2007| 
have discussed the utility of the Fourier method (Gray 
I1976P in the case of early type stars. Each of these meth- 
ods have their advantages and disadvantages, with for ex- 
ample, the Fourier method being able to deconvolve the 
rotational broadening from other broadening mechanisms 
such as macroturbulence, although it requires high quality 
spectra. Profile fitting to obtain the rotational velocity is 
possible for lower quality spectra, but relies on an appro- 
priate description of the intrinsic line profile (from a stellar 
atmosphere code), and the broadening agents affecting the 
line profile. Measuring the FWHM and using the cross cor- 
relating method is relatively straight forward but does not 
take into account the possible extra non-rotational broad- 
enings (such as Stark broadening of the H and He lines, 
macroturbulence and microturbulence) thereby leading in 
some cases to only upper limits for the projected rotational 
velocity estimate. 

Given the quality of our spectra we decided to use the 
profile fitting methodology discussed in Paper III for all 
the non-supergiant objects as the quality of our data is 
not high enough to make the Fourier method viable for 
the entire sample. Briefly, a model spectrum, at the closest 
TLUSTY grid point to the derived atmospheric parameters, 
was selected and the equivalent width of an appropriate 
line was then scaled to the same strength as that observed. 
Typically this scaling was less than 10% and, except for 
the narrowest lined spectra, has little effect on the derived 
projected rotational velocity. The model profile was first 
convolved with a Gaussian profile to take into account the 
spectral resolution of our data; then, the resultant profile 
was convolved with a rotational broadening function for a 
range of rotational velocities. A chi-squared minimisation 
test was used to estimate the projected rotational velocity. 
The atmospheric parameters were then redetermined with 
this new estimate and the process repeated. Generally no 
more than two iterations were required. 

Hydrogen lines were observed in all our spectra, but 
their profiles are dominated by Stark broadening. Although 
this also affects the diffuse He i lines, at moderate projected 
rotational velocities, the rotational broadening begins to 
dominate, and hence they are suitable for estimating this 
quantity. The Hei 4026A line is observed across all our 
spectral types and additionally is well separated from other 
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absorption features. Therefore this hne was employed for 
the profile fitting methodology. Note however, that at low 
projected rotational velocities the profile of this He i line is 
no longer sensitive to the rotational broadening. For these 
cases, we have instead used metal lines if available (typi- 
cally either the Mgii 448lA line or the Sim 4552A line). 
These projected rotational velocity estimates are listed in 
Tables |3]l6l apart from those of the supergiant targets dis- 
cussed below. 

A small fraction of the objects in our sample have 
uncertain spectral types that lie in the range B 1-3 and, 
as discussed above, the He line spectrum is not capable 
of constraining their temperatures. In order to derive the 
projected rotational velocity of these objects, approximate 
models have been adopted. As these stars frequently have 
large rotational velocities (leading to the uncertain spectral 
classification), the broadening of the observed lines will be 
rotationally dominated and hence assuming an inappropri- 
ate model will not significantly affect the derived rotational 
velocity. Atmospheric parameters are not given in Tables 
[S]for these objects. 

We have also estimated rotational velocities using the 
Fourier transform methodology where possible. A detailed 
discussion of the applicability of this methodology in the 
case of early type stars is given by Simon-Diaz & Herrero 
(|2007p and will not be repeated here. In Fig. [1] we compare 
the projected rotational velocity derived by the profile fit- 
ting methodology and the Fourier transform methodology 
for the non-supergiant targets, where the spectra were of 
sufficient quality for the application of the Fourier method. 
We find in almost all cases good agreement (within 10%), 
indicating that the choice of methodology is unimportant 
for these targets and validating our adopted projected ro- 
tational velocity estimates. 



illustrated by Ryans et al (PUO^ . Dufton et al. (|2006bp and 
Simon-Di'az & Herrero (|2007p . there is an important non- 
rotational broadening mechanism affecting the profiles of 
early type supergiants (usually termed macroturbulence) . 
Since the line fitting methodology does not allow for this ex- 
tra source of broadening, the Fourier transform technique 
has been used to derive the projected rotational velocity 
of objects classed as supergiants in the sample. In Fig. [5] 
we compare the derived rotational velocities from the two 
methods for these objects. Due to the quality of the spec- 
tra, in terms of signal to noise ratio, in many cases it was 
only possible to derive upper limits to the projected rota- 
tional velocity from the Fourier method (see Simon-Diaz & 
Herrero); these are indicated by downward pointing arrows 
in Fig. [1 
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Fig. 1. Comparison of the projected rotational velocity de- 
rived by the profile fitting method and the Fourier trans- 
form method for the non-supergiant objects in the sam- 
ple. The dotted line indicates a one-to-one correlation. 
The dashed lines indicate a 10% or lOkms^^ uncertainty, 
whichever is the larger. 



As already proposed by several authors (Slettebak ll956l 
Conti & Ebbets fWTl Howarth et al. [TM7|) . and recently 



Fig. 2. Comparison of the projected rotational velocity 
(wsini) derived by the profile fitting method and the 
Fourier transform method for the supergiant objects in 
the sample. The dotted line indicates a one-to-one corre- 
lation. The dashed lines indicate a lOkms"^ uncertainty. 
Downward pointing arrows indicate those objects where it 
was only possible to derive upper limits to the v sin i from 
the Fourier method. 

From Fig. [5] it is apparent that in many cases the 
projected rotational velocities derived from the Fourier 
Transform method are significantly lower than those from 
the profile fitting method and these are likely to be objects 
with significant macroturbulent broadening. Hence for the 
supergiant targets, the projected rotational velocities listed 
in Tables are those deduced from the Fourier Transform 
methodology. We note that for the determination of the at- 
mospheric parameters of the supergiant objects in the sam- 
ple, the estimated rotational velocity determined from the 
profile fitting method was used rather than that from the 
Fourier Transform method. This is valid as it is the entire 
broadening of the line which is important when determining 
the parameters and not simply the rotational component. 

The uncertainties in effective temperature and surface 
gravity have little effect upon our estimated projected rota- 
tional velocities. For example, NGC 346-083 was analysed 
with a model profile at a temperature of 27 500 K and a 
surface gravity of 4.00 dex and a projected rotational ve- 
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locity of 207kms~^ was estimated. Adopting a higher tem- 
perature profile of 30 000 K and higher surface gravity of 
4.25 dex (as the surface gravity and effective temperature 
estimates are directly correlated) leads to a w sin i estimate 
of only 4kms~^ greater. Similarly adopting a lower tem- 
perature and gravity (25000 K, 3.75 dex) results in a value 
14kms~-'^ greater. Including the effects of normalisation, 
scaling and choice of model profile we would not expect the 
errors in our rotational velocity estimates to exceed 10%. 
This is confirmed by the comparison shown in Fig. [1] 
where different lines were used for the Fourier Transform 
and the profile fitting methods. We note that the Fourier 
Transform method is additionally model independent. 

Rapidly rotating objects (close to critical velocity) can 
not be considered spherical and equatorial gravity darken- 
ing (the cooler and fainter equatorial regions of an oblate 
star contribute less flux to the observed spectrum) can sig- 
nificantly alter the estimated atmospheric parameters and 
projected rotational velocities (von Zeipel|T924|). Townsend 
et al. (|2004p have shown that this effect leads to an un- 
derestimation of projected rotational velocities for close to 
critical rotational velocities and that the maximum ratio of 
projected rotational velocity to critical velocity {vsini/vc) 
that can be inferred from line widths is approximately 0.8. 
We have derived critical rotational velocities for our sample 
and find that less than 5% of our sample have v sini/uc 
ratios greater than 0.5 with the maximum value of this 
ratio being approximately 0.7 (NGC 2004-113). If we con- 
sider that this object is rotating at near to critical velocity, 
the projected rotational velocity would be underestimated 
by less than 15% due to the effects of gravity darkening 
(Fig. 1, Townsend et al.). Hence we believe that gravity 
darkening effects should be small compared to our observa- 
tional uncertainties for the majority of our sample of stars 
and in the most extreme cases of the same order. 

3.3.1. Comparison with other work 

Martayan et al. (|2006|, I2007P have observed a large number 
of Be and B-type stars towards NGC 2004 and NGC 330 
also with the FLAMES instrument, but with lower spec- 
tral resolution than we employed. In order to make com- 
parison with the results of Martayan et al. the samples 
were constrained to the magnitude range where there is 
a good overlap between them. For NGC 330 this is the vi- 
sual magnitude range of 15.0-16.5 with a comparison of 
the two samples being shown in Fig. [S] The cumu- 
lative probability distributions of the B-type stars closely 
agree. There is some suggestion that the FLAMES sample 
contains a larger fraction of stars with very low projected 
rotational velocities than does the that of Martayan et al. 
Given the higher spectral resolution and signal to noise ra- 
tio of our FLAMES data, this is not unexpected. There 
are 13 objects in common between our FLAMES survey of 
NGC 330 and that of Martayan et al. and in Fig. 2] we com- 
pare the projected rotational velocities for these stars. In 
general there is reasonable agreement between the two sam- 
ples, although the discrepancy at low projected rotational 
velocities is again evident. 

The distributions of the Be stars shown in Fig. [3] are in 
clear contrast, with Martayan et al. (2007) finding system- 
atically higher projected rotational velocities. The reason 
for such a large discrepancy is unclear. For the estimation of 
parameters of their Be sample, Martayan et al. removed the 




V sin i / km s 



Fig. 3. Comparison of our estimated rotational velocities 
(FLAMES) with those of Martayan et al. ()Wf[ Mar07) for 
the B and Be type stars towards NGC 330. This clearly 
shows the discrepancy between the Be analysis pre- 
sented here and that of Martayan et al. Only objects 
in the visual magnitude range 15-16.5 have been included. 




V sin i - FLAMES (l<m/s) 

Fig. 4. Comparison of our estimated rotational velocities 
(FLAMES) with those of Martayan et al. (2007, Mar06). 
The solid line indicates a one-to-one correlation. The 
dashed lines indicate our 10% uncertainty down to a limit 
of 10kms~\ 



core of the lines where Be emission was observed and hence 
the projected rotational velocities are estimated from the 
wings of the line. In Fig. [5] we present a comparison between 
a model spectrum of a H7 line with no rotational broaden- 
ing and that with a rotational broadening of 300kms~^. It 
is clear that the wings are not as sensitive as the line core 
to changes in the rotational broadening. As such removing 
the line core will increase the uncertainty of the rotational 
velocity estimate, especially for low signal to noise data. 
However, by normally using only one line to determine the 
projected rotational velocity, our estimates may be biased 
towards smaller rotational velocities due to saturation ef- 
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fects (Fremat et al. I2005P compared to those of Martayan 
et al. who have used several lines, although, this should 
only be an issue at high rotational velocity. As discussed 
above gravity darkening should be negligible for the ob- 
jects in our sample. Indeed, for ten Be-type objects we 
have also derived rotational velocities from the He I line at 
4387A and generally find good agreement with the values 
from the 4026A line. Additionally in estimating the pro- 
jected rotational velocities, models with a surface gravity 
of 4.0 dex have been adopted. At this high gravity the Stark 
broadening will be relatively large and hence less rotational 
broadening is required to fit the profile compared to assum- 
ing a lower gravity, see Table El As such our Be star rota- 
tional velocities may be biased towards lower values, but 
this is only an issue at low projected rotational velocities 
and in these cases we have endeavored to use metal lines 
where possible in order to minimize this bias. 



stars are typically the fastest rotators in the sample, this 
could bias our velocity distribution towards lower rotational 
velocities. To investigate this we ran the target selection 
software, configure, again for the NGC330 field, with- 
out explicitly excluding any objects from the input cat- 
alogue. From this test we found seven of the previously 
excluded stars were included in the resulting fibre config- 
uration (which of course differs from our actual observed 
field by more than seven stars, due to the different posi- 
tions of the allocated fibres). If we consider that the net 
difference was to replace seven normal B-type stars with 
seven Be targets, the mean rotational velocity of the SMC 
stars would only increase by ~3kms~^. Given that some 
of the excluded objects are not actually confirmed Be stars 
this is a conservative estimate, and we believe that any mi- 
nor selection effect that may have been introduced does not 
unduly bias our conclusions. 



1 .0 



2. 0.8 



o 
z 

0.4 




4320 4330 4340 4350 4360 
Wavelength (angstroms) 

Fig. 5. Comparison of a theoretical II7 line with no rota- 
tional broadening (solid line) with the same line rotation- 
ally broadened by 300kms~^ (dashed line) showing that 
the wings of the hydrogen line are not sensitive to 
rotational broadening. 

We have no NGC 2004 objects in common with the sam- 
ple of Martayan et al. (i2006|) and the magnitude range of 
overlap is much smaller than that of the NGC 330 samples 
although similar trends are seen in the rotational veloc- 
ity distributions. The B-type stars are in excellent agree- 
ment and again there is some evidence that the Be stars of 
Martayan et al. have higher rotational velocities. Although 
the discrepancy between the Be star analyses cannot be 
fully resolved, this does not invalidate the metallicity de- 
pendence discussed in Sect. 15.21 Indeed, since the SMC 
has a higher fraction of Be stars than the LMC, using 
the results of Martayan et al. would enhance the differ- 
ence between the rotational velocity distributions of the 
two Magellanic Clouds. 

3.3.2. Selection effects 

As discussed in Paper II approximately 30 (potential) Be 
stars were removed from the list of possible targets for 
the FLAMES pointing towards NGC 330. Given that Be 



3.4. 0-type stars 

Atmospheric parameters and projected rotational velocities 
of the majority of the 0-type objects observed Paper II 
have been derived by Mokiem et al. f20061 '2007) and these 
have been directly inserted into Tables [3][6l We note that 
for those stars where our analysis overlaps with that of 
Mokiem et al. we find good agreement between the two 
sets of atmospheric parameters and rotational velocities. 

3.5. Binary objects 

After correcting for radial velocity shifts (see Sect. 12. the 
single lined spectroscopic binaries were analysed using the 
methodology outlined above. However, this was not possi- 
ble for the double lined spectroscopic binaries, due to both 
continuum contamination and the double line nature of the 
profiles. In Tables we indicate those objects for which 
we see evidence of a secondary object, although for small 
asymmetries and velocity shifts in the line profile, variabil- 
ity may also be an explanation. For all cases, the individual 
exposures were examined to find those in which the line pro- 
files from the primary and secondary objects had the max- 
imum separation. The profile fitting methodology was then 
used to estimate the rotational velocity of the primary, and 
where possible, the secondary object. We do not attempt 
to derive atmospheric parameters for either object. Metal 
lines were used to estimate the rotational velocities in order 
to minimise the uncertainty from the choice of model pro- 
file. An example of the fit for the double lined spectroscopic 
binary system NGC 346-020 is shown in Fig. [51 where the 
Sim lines from each object are partially blended. Where 
it has been possible to estimate the projected rotational 
velocity of the secondary object in a binary system this 
is listed in Tables [3]|6] with the letter 'B' appended to the 
system's identifier. 

3.6. Be objects 

A significant fraction of the stellar sample is made up of 
Be type objects, with the majority having spectral types in 
the range B 1-2. In Tables we do not give atmospheric 
parameters for these objects due to problems associated 
with continuum contamination from the circumstellar disc 
and emission lines. Additionally as many of the objects do 
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Table 7. Dependence of the derived usini value upon surface gravity for a sample of Be objects. For large rotational 
velocities, the Sim lines could not be identified. 



Line log g usini (kms ^ 

(dex) NGC 346-091 NGC 346-095 NGC 346-096 

He I 4026A lOO 74 227 346 

Hei 4026A 3.50 120 250 343 

Siiii 4552A 4.00 49 

Siiii 4552A 3.50 47 




o 



0.94 - 

0.92 \ ^ 

4545 4550 4555 4560 

Wavelength (Angstroms) 

Fig. 6. Fit of the Si ill line at 4552A for the double Hned 
spectroscopic binary system NGC 346-020. The model spec- 
tra for the primary and secondary objects have been fitted 
with a vsini value of 34kms~^ and 66kms~^ respectively. 



not have precise spectral types it was not possible to assign 
an effective temperature based on spectral type. We have 
however used the metal line and helium spectra to estimate 
the projected rotational broadening. 

In Table[7]we show how the estimated projected rota- 
tional velocity depends on the adopted surface gravity for 
three Be type stars with low, medium and high velocities. 
It is evident that at large rotational velocities the He profile 
is dominated by the rotational broadening and is no longer 
sensitive to the adopted profile. Additionally the projected 
rotational velocity estimated from the metal line spectra 
does not appear to be dependent on the adopted model 
profile and we have adopted these estimates in Tables [3]l6] 
where possible. 

3.7. Luminosities 

Luminosities were calculated for each object in the sample 
using the same methodology as discussed in Hunter et al. 
(|2007p . Constant values of reddening towards each cluster 
were assumed, with E{B — V) values of 0.09 (Massey et al. 
fT995l) . 0.06 (Gonzalez & Wallerst ein gM Sl), 0.13 (Massey 
et al.) and 0.09 (Sagar & Richtler flQQip being adopted for 
NGC 346, NGC 330, Nil and NGC 2004 respectively. We 
adopt the SM C redd ening law of = 2.72E{B - V) from 
Bouchet et al. ([T985)) for the SMC targets while we use the 
standard Galactic law for the LMC targets. We note that 



assuming the standard Galactic law would change the de- 
rived SMC luminosities by less than 0.04 dex. Bolometric 
corrections have been adopted from Vacca et al. ( I1996P 
for objects hotter than 28,000 K and Balona (1994]) 
for cooler objects, together with distance modulii of 
18.91 dex (Hilditch et al. 2005) and 18.56 dex (Gieren et al. 
[2005) for the SMC and LMC. The luminosities are given in 
Tables and typically have an uncertainty of 0.1 dex. 

3.8. Masses 

In Fig. [71 Hertzsprung-Russell diagrams are presented for 
the sample of stars towards each cluster. Non-rotating evo- 
lutionary tracks are also shown and from these the evo- 
lutionary masses of each object have been estimated and 
these are listed in Tables ISHll Isochrones are not included as 
our sample is likely to be dominated by field stars and there- 
fore the populations are unlikely to be coeval, see Sect. 12.31 

The dashed line in Fig. [7] indicates the magnitude limit 
for our observations in each region (see Paper II) and it is 
clear that we sample the zero-age main-sequence (ZAMS) 
at differing masses in each clusters. For example we tend to 
sample more evolved objects in NGC 2004 due to our target 
slection, but do sample stars on the ZAMS below 9, 12 and 
15 M0 in NGC 330, NGC 346 and Nil respectively. 

The uncertainty in the masses for those objects for 
which we have estimated effective temperatures is typically 
10-15%, which arises principally from the uncertainty in 
the luminosity. In Tables we do not give atmospheric 
parameters for Be stars, double lined spectroscopic bina- 
ries and stars with uncertain spectral types. However, for 
the estimation of their masses, we assume temperatures 
based on their spectral type (or their mean spectral type 
where a spectral type range is given). The uncertainty in 
the mass of these objects is typically 30%. In addition to 
these evolutionary masses it is possible to estimate spectro- 
scopic masses, but an uncertainty of 0.25 dex in the surface 
gravity estimate corresponds to an uncertainty of over 70% 
in the spectroscopic mass and hence such estimates are of 
limited value. 

3.8.1. Mass discrepancy 

Although individual spectroscopic masses are highly un- 
certain, our sample is large enough to examine systematic 
differences between evolutionary and spectroscopic masses. 
For example, it is well known that the evolutionary masses 
derived for OB-type supergiants typically are greater than 
the spectroscopic masses (see for example, Herrero et al., 
1992; Trundle & Lennon, 2005; Lennon et al., In 
Table [8] the average mass discrepancies are summarized. 
While the discrepancy for evolved stars has previously been 
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Fig. 7. Hertzsprung-Russell diagrams for our sample of objects towards Nil, NGC 2004, NGC346 and NGC330. Be 
objects are plotted as squares and objects showing evidence of binarity are plotted as open symbols. Non-rotating 
evolutionary tracks have been obtained from Meynet et al. p994p and Schaerer et al. (|1993p for the LMC objects and 
from Meynet et al. and Charbonnel et al. p993p for the SMC objects. The dashed line indicates the magnitude limit for 
our observations in each region. It should be noted that many of the objects in each sample lie in overlapping positions 
in these figures. 



observed, that for the main-sequence objects is unexpected. 
Although the scatter is large, given the number of objects, 
there may be a systematic difference. If it is assumed that 
no mass discrepancy should exist for main-sequence objects 
since this evolutionary stage is thought to be well under- 
stood, that observed could be due to binarity, errors in the 
distance estimate, bolometric corrections, reddening values 
or even the derived surface gravities. However, given that 
the magnitude of the effect is different in the SMC com- 
pared to the LMC, it is unlikely to be due to the derived 
gravities, or bolometric corrections as the same methodolo- 
gies were used for each. Additionally the mass discrepan- 
cies for the binary systems are almost identical to those 
in Table [S] indicating that undetected binary systems are 
unlikely to be the cause of the discrepancies. Since the red- 
dening values are small, the discrepancy may come from 
the distance estimates. 

Using the evolutionary masses and atmospheric param- 
eters listed in Tables [3][6l we have derived luminosity es- 



Table 8. Mean mass-discrepancies for the sample of sin- 
gle B-type stars. Objects are classified by surface 
gravity with supergiants having gravities less than 
3.2 dex (see Sect. 15. ip . giants between 3.2 and 
3.7 dex and dwarfs greater than 3.7 dex. Errors are 
the standard deviations and the value in brackets indicates 
the number of objects. 



log (IVIcvolutionary /-^^spectroscopic ) 

LMC SMC 

Supergiants 0.17±0.17 (22) 0.25±0.19 (13) 
Giants 0.06±0.14 (43) 0.18±0.14 (39) 

Dwarfs -0.14±0.22 (51) -0.05±0.17 (60) 



timates. It should be noted that these are different from 
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the luminosities listed in Tables which assumed a fixed 
distance to each Cloud. Adopting the same reddenings and 
bolometric corrections as discussed in Sect. 13. 7[ it is then 
possible to derive distance estimates towards each star. For 
objects with logg >3.7dex the mean distance moduli to 
the LMC and SMC are 18.16±0.77 and 18.79±0.55 from 
51 and 60 objects respectively, where the imcertainties are 
the la standard deviations, compared with our adopted 
values of 18.56 and 18.91. Although the mean value for the 
LMC is 0.4 dex smaller than that adopted in Sect. 13.71 it 
lies within the range of those reported for the LMC in the 
recent literature. For example, Udalski (1998) estimates a 
distance modulus of 18.13 for the LMC. However, it should 
also be noted that Udalski similarly estimates a distance 
modulus of 18.63 dex for the SMC which is 0.17 dex closer 
than our derived value of 18.80 dex and 0.28 dex closer 
than our adopted value of 18.91 dex. Excluding the 
most massive stars (greater than 25 q) increases the 
mass discrepancies although within the uncertain- 
ties this is not significant. To maintain consistency with 
published work from the FLAMES project the distance esti- 
mates discussed in Sect. ISTfl have been adopted throughout 
the subsequent discussion. 

If corrections are applied to remove the mass discrep- 
ancy of the dwarf stars it further enhances the mass dis- 
crepancy of the evolved objects. However, such corrections 
result in similar mass discrepancies for the LMC and SMC 
evolved samples. This may indicate that the post-main- 
sequence evolution of massive stars is not fully understood 
and as discussed above, such mass discrepancies have pre- 
viously been observed. In Sect. 15.11 several inconsistencies 
between the observations and the theoretical predictions 
of evolved objects are highlighted. In particular, we argue 
for extension of the hydrogen burning phase of the evo- 
lutionary tracks which would result in lower evolutionary 
masses being derived, hence reducing the mass discrepancy 
for evolved stars. 

4. Populations within the sample 

Atmospheric parameters and rotational velocities 
have been derived for a large sample of objects 
which contains several groups with different prop- 
erties. In this section we compare these various sub- 
sets using as a control sample single, core hydrogen 
burning stars of less than 25 Mq. In Sect. [5] we dis- 
cuss how the subsets can be utlised to answer the 
three questions listed in Sect. [TJ 

4.1. High mass stars 

Mass-loss rates generally increase with stellar mass (Maeder 
& Meynet I200ip and hence more massive objects should 
lose more angular momentum and rotate slower than less 
massive objects for a given initial rotational velocity. The 
theoretical evolutionary tracks of Maeder & Meynet show 
that mass-loss effects become significant for masses greater 
than 25 Mq. In TablelHlwe compare the mean projected ro- 
tational velocity for stars with masses greater than 25 Mq 
and those with masses of < 25 Mq . There is indeed some ev- 
idence that the more massive stars rotate slower, although 
it should be stressed that the number of objects with masses 
greater than 25 Mq is only 7% of the sample and the 
trend may not be statistically significant. 



Table 9. Mean projected rotational velocities for single 
core-hydrogen buring stars with masses greater than 
25 Mq and masses of < 25 Mq in each region. Values in 
brackets indicate the number of stars in each sample. 



Region Mean vs'mi (kms '") 
M<25M0 M>25Mq 



NGC 346 


150 


(65) 


113 


(5) 


NGC 330 


154 


(77) 


73 


( 1) 


Nil 


149 


(34) 


126 


(12) 


NGC 2004 


121 


(60) 




(0) 



4.2. Binarity 

Binarity effects can significantly alter the evolution of mas- 
sive stars. For example, mass-transfer can lead to spin-up 
of the accreting star, while tidal interaction in close binary 
systems can lead to reduced rotational velocities through 
partial or full synchronisation of the rotational velocity with 
the orbital velocity (see for example Huang & Gies 120061 
Abt et al.[inn3Zahn[T177|). Indeed both Zahn ()TMi|) and 
Hunter et al. (|2007p find that massive stars in binary sys- 
tems typically have smaller amounts of core processed ma- 
terial in their photosphere compared to single stars. This 
is consistent with tidal interaction reducing the rotational 
velocities and hence suppressing the rotational mixing of 
the products of nucleosynthesis into the photosphere. 

Table 10. Mean projected rotational velocities of 
the core hydrogen burning single and binary stars 
less massive than 25 Mq in each region. The binary 
population includes both single and double lined 
spectroscopic binaries. Values in brackets indicate the 
number of stars in each sample. 



Region Mean vs'mi (kms ) 

Single stars Binary stars 



NGC 346 


147 


(67) 


103 


(25) 


NGC 330 


146 


(83) 


121 


(8) 


Nil 


143 


(35) 


92 


(30) 


NGC 2004 


120 


(62) 


99 


(28) 



In Table [TOl the mean projected rotational veloc- 
ities of single and binary stars are compared. It can 
be seen that binaries tend to have lower rotational 
velocities than single stars. This may indicate that 
tidal forces do play an important part in reducing 
the rotational velocities of massive stars. However, 
it is important to consider selection effects. Radial velocity 
variations (the evidence of binarity) are detected by cross- 
correlating individual exposures of the same wavelength re- 
gion observed at different epochs to look for radial velocity 
shifts (Sect. [^TT|) . Such shifts are relatively simple to detect 
to a few km s^^ for stars with projected rotational velocities 
up to 100 km s^^ as the absorption lines are typically sharp 
and well defined. For stars with broad absorption lines and 
hence large projected rotational velocities, it is difficult to 
define the position of the core of the line and small radial ve- 
locity shifts will remain undetected. Additionally in binary 
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Table 11. Mean projected rotational velocities of the sin- 
gle supergiant (surface gravities of less than 3.2 dex) 
and core hydrogen burning stars with masses <25M0. 
Values in brackets indicate the number of stars in each sam- 
ple. Note that as upper limits are included in the supergiant 
sample, the means given here should be treated as upper 
limits to this class of object. 



Region 


Mean vsini (kms~^) 




Supergiants Core hydrogen 




burning stars 



NGC 346 


72 


( 1) 


147 


(67) 


NGC 330 


38 


(12) 


146 


(83) 


Nil 


60 


(11) 


143 


(35) 


NGC 2004 


60 


(8) 


119 


(62) 



systems with small orbital periods, the separation between 
the objects would be smaller, the velocity shifts greater 
and tidal interaction would be more important than for 
longer period systems. As such the rotational velocity dis- 
tributions of binary objects could be biased towards stars 
having low rotational velocities. Clearly further studies of 
the binary populations in young clusters are required. 

4.3. Supergiants 

In Table 1111 the mean projected rotational veloc- 
ities of the supergiant and core hydrogen buring 
objects are compared with supergiants being con- 
sidered as objects with surface gravities less than 
3.2 dex (see Sect. [STT]) . Assuming that the inclina- 
tion axes of our targets are randomly orientated the 
mean rotational velocities would be a factor of 4/7r 
larger than the values quoted in Table [TTl It is clear 
that supergiants have much lower rotational veloc- 
ities than core hydrogen burning objects which is 
simply a consequence of their evolutionary status. 

4.4. Be stars 

It has recently been reported that Be stars typically have 
larger rotational velocities than B-type stars (Martayan et 
aL [M7Sl[Myf|) . A significant number (17%) of the stars in 
the FLAMES survey have been classified as Be type objects 
and hence it is possible to compare the rotational velocities 
of normal B-type stars with those showing Be characteris- 
tics. In Table[T2]we present the average projected rotational 
velocities of the B- and Be-type stars in each region. 

It is apparent that on average Be-type stars have higher 
rotational velocities than normal B-type stars in agreement 
with the finding of Martayan et al. (|2006[ [2007|) . see also 
Fig. [3] Zorec et al. (|2005p have shown that the Be phe- 
nomena normally occurs early in the core hydrogen burn- 
ing phase. They also show that for the most massive stars 
the Be phenomena disappears during their main-sequence 
lifetimes and postulate that this is due to mass-loss causing 
a star to spin down. 

Table[T2]additionally shows that we have observed many 
more Be-type stars in the SMC than the LMC, again in 
agreement with Martayan et al. 120071 If B-type stars be- 
come Be-type stars because they reach the main-sequence 



Table 12. Comparison of the mean projected rota- 
tional velocities of the single core hydrogen burning 
B- and Be-type stars with masses <25Mq in each 
region. Values in brackets indicate the number of stars in 
each sample. 



Region Mean wsini (kms 'y 

B stars Be stars 



NGC 346 


115 


(45) 


213 


(22) 


NGC 330 


128 


(64) 


207 


(19) 


Nil 


134 


(31) 


215 


(4) 


NGC 2004 


105 


(51) 


190 


(11) 



with sufficiently high angular momentum, this may indi- 
cate that at lower metallicity, massive stars tend to rotate 
faster and suffer a smaller amount of spin down. 

5. Discussion 

5.1. Constraints on evolutionary history 

The surface gravity of an object is a good indicator of its 
evolutionary status, with the gravity decreasing as an ob- 
ject evolves from a dwarf, to a giant and then into a su- 
pergiant star. In Fig. [5] we plot the projected rotational 
velocity against surface gravity. 

The distributions of rotational velocities shown in Fig. [5] 
appear bi-modal in the sense that for surface gravities of 
<3.2dex the projected rotational velocities are small, while 
at higher values the distribution is much broader. We in- 
terpret this limit of 3.2 dex to be the end of the 
core hydrogen burning phase, since after this point the 
stars rapidly expand to become slowly rotating, core helium 
burning objects, as represented by the late B-type objects 
in our sample with surface gravities of ~2.0dex. Hence we 
consider objects with surface gravities lower than 
3.2 dex to be supergiants. We note that the predicted 
position of the end of core hydrogen burning shown by the 
evolutionary tracks occurs at higher surface gravities than 
that implied by our observations and suggests that the pre- 
dicted core hydrogen burning lifetime needs to be extended 
to lower surface gravities. It is however important to stress 
a caveat. If the mass-discrepancy discussed in Sect. 13.8.11 
is not due to an error in the adopted distance moduli, but 
instead due to a systematic error in our estimated surface 
gravities then our targets would have surface gravities that 
are systematically '^0.1 dex higher. 

A large population of the giants lie in the post-main- 
sequence gap when compared with the evolutionary tracks 
shown in Fig. [71 Due to their small evolutionary timescale, 
such objects would not normally be expected to be ob- 
served. A similar over abundance of post- main-sequence ob- 
jects in the Galactic clusters NGC 3293 and NGC 4755 was 
observed in Paper III. Extension of the end of the main- 
sequence (for example, by increasing the efficiency of over- 
shooting) would allow the evolutionary tracks to encompass 
this giant phase during the core hydrogen burning phase. 
This would then be consistent with the presence of such 
objects in our samples. 

The theoretical predictions of Maeder & Meynet l|200ip 
indicate that rotational velocities remain relatively con- 
stant during core hydrogen burning, decreasing by approx- 
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Fig. 8. Plots of projected rotational velocity (usini) against the surface gravity for stars towards Nil, NGC2004, 
NGC 346 and NGC 330 showing the bi-modal distribution of rotational velocities. Be objects are assumed to be 
main-sequence objects with a surface gravity of 4.00 dex and are plotted as diamonds. Single lined spectroscopic binaries 
are indicated by square symbols. Downward arrows indicate upper limits to the v sini. The stellar evolutionary tracks of 
Maeder & Meynet (|200ip for a 15 Mq star, with an initial rotational velocity of 300 km s"""^ (scaled to take into account 
random angles of inclination) at SMC (dashed line) and Galactic (solid line) metallicity are shown. 



imately 40kms^^ from the zero-age main-sequence to the 
end of hydrogen burning. From Fig. [5] there may be some 
evidence that the objects with the largest surface gravities 
cover the widest range of rotational velocities and indeed 
objects classified as giants have rotational veloci- 
ties that are lower than those of dwarfs by on aver- 
age ~50kms~^ for our entire sample. This appears to be 
consistent with giants being in the core hydrogen burning 
phase, again implying that extension of the evolutionary 
tracks may be necessary. 

The evolutionary tracks of Maeder & Meynet (I200ip 
which incorporate rotation show that for a 15 Mq star it 
takes ~13.6Myr to reach the end of the hydrogen burning 
main-sequence, while the blue supergiant phase has a life- 
time of less than 0.1 Myr, i.e. the blue supergiant phase 
is less than 1% of the core hydrogen burning life- 
time. The large number of supergiants observed in 
the FLAMES survey is unexpected with the ratio of 
supergiants to core hydrogen burning objects being 
~11%. Significant populations of objects in the post-main- 
sequence gap have previously been seen (see, for example, 
Caloi et al. ll993l Fitzpatrick & Garmanv [T990)) . However, it 
should be noted that the FLAMES survey preferentially se- 



lected the brightest objects towards each cluster, and hence 
there may be some bias towards selecting supergiants. 

The evolutionary tracks shown in Fig. [5] imply that the 
rotational velocity slows down relatively smoothly as the 
star evolves to low surface gravities after the core hydrogen 
burning phase, although, as discussed above, this phase 
is almost instantaneous compared with the core hydrogen 
burning lifetime. However, observationally this trend is not 
seen; objects with surface gravities <3.2dex have relatively 
constant rotational velocities, which are smaller than pre- 
dicted. A similar trend has previously been observed by 
Dufton et al. (|2006bp . However, as discussed by Hunter et 
al. (|2007p and Trundle et al. (2007) many of the supergiant 
objects are highly enriched in nitrogen and show evidence 
of binarity, suggesting that a mass-transfer event may have 
occurred. As such, comparison with single star evolution- 
ary tracks may not be valid. Indeed, rejuvenation through a 
mass-transfer event is one method of populating this short 
evolutionary stage and may reconcile the over-abundance 
of supergiants in our samples with the predicted lifetime of 
such an evolutionary stage from single star models. 

An alternative explanation is the possibility of blue- 
loops, where the supergiant goes through a red-supergiant 
phase and then returns to higher effective temperatures. 
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Fig. 9. Rotational velocity distributions of the single core-hydrogen burning LMC and SMC stars with masses 
<25Mq, upper and lower panels respectively. The LMC distribution is fitted with a intrinsic rotational velocity 
modelled by a Gaussian with a peak at 100 kms^"'^ and a ^ half width, of 150kms~^, while the SMC distribution is fitted 
with a Gaussian having a peak at 175kms~^ and a half- width of 150knis^^. 



Such a process would lead to the observed nitrogen en- 
hancements and rotational velocities and an increased life- 
time in the blue supergiant phase although the blue-loop 
phase is not well constrained in the evolutionary models 
(G. Meynet and N. Langer, private communications). We 
therefore beheve that the supergiant population of objects 
shown in Fig. [8] are not the result of direct single star evolu- 
tion from the core hydrogen burning phase and either 
binarity or blue loops may be responsible for this 
supergiant population. 

5.2. Rotation and metallicity 

As discussed in Sect. [1] stars at low metallicity are pre- 
dicted to rotate faster as they lose less angular momentum 
through their stellar winds and are generally more compact. 
Our sample of objects covers the metallicity regimes of the 
LMC and SMC and hence it is possible to look for metal- 
licity effects. To increase our sample sizes we have com- 
bined together those for Nil and NGC 2004 in the LMC 



and NGC 346 and NGC 330 in the SMC sample. Only sin- 
gle core hydrogen burning objects with masses of 
<25Mq have been considered and this corresponds 
to 48% and 70% of our LMC and SMC samples 
respectively. The upper limit in mass was set to 
minimise the possible effects due to mass-loss as 
discussed in Sect. 14. H whilst the supergiants have 
been excluded as their rotational velocities will not 
be representative of their main-sequence progeni- 
tors. Our samples generally have spectral types be- 
tween BO and B3, and have initial masses between 7 
and 25 Mq. The spectral type distributions of both 
samples are similar and there is no reason to expect 
any trend of rotational velocity with spectral type 
over this small range. 

In Fig. [9l the distribution of projected rotational veloc- 
ities in the LMC and SMC arc shown together with their 
cumulative probability functions. The number of objects we 
have observed should be sufficient to allow us to estimate 
the intrinsic rotational velocity distribution. Following the 
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procedures discussed in Paper III, we have attempted to fit 
these distributions with a Gaussian function convolved with 
a function to take into account random angles of inclination. 
A chi-squared test has been performed to obtain the best 
fit to the cumulative probability functions. We chose not to 
apply the chi-squared test to the fit to the histograms, as 
this may be dependent on the choice of binsize. 

For the LMC sample, an intrinsic velocity distri- 
bution with a peak at 100 km s"^ and a ^ half width 
of 150 km s^^ (which corresponds to a to a full-width- 
half-maximum, FWHM, of 250 km s^^) is found, with 
corresponding values for the SMC sample being 
175 km s^^ and 150 km s^^. Although a Gaussian dis- 
tribution has been assumed, as discussed in Paper III and 
in Mokiem et al. (|2006p any function with a similar pro- 
file would also provide a satisfactory fit. For example as 
the Gaussian profiles are quite broad, a rectangular pro- 
file centred on the Gaussian peaks, with similar widths, 
would also give a reasonable fit to the observed distribu- 
tions. From these results it appears that the LMC objects 
typically rotate significantly more slowly than those of the 
SMC. Indeed, the Student's t-test gives a 93% probabihty 
that the means of the two distributions are significantly 
different while the K-S statistical test gives the probabil- 
ity that the distributions are drawn from the same parent 
population to be less than 10%. 

To further investigate metallicity effects the Magellanic 
Cloud samples can be compared with Galactic results. 
However before this is possible, it is necessary to ensure 
that the samples are comparable. For example, the Galactic 
objects observed in the FLAMES project and discussed in 
Paper III mostly lie within two cluster radii of their respec- 
tive clusters. From the discussion in Sect. I2.3l it is apparent 
that the objects that we observed in the Magellanic Clouds 
lie much further from the centre of their respective associ- 
ations. Our sample is therefore likely to be dominated by 
field or unbound objects and it may be more appropriate to 
compare the rotational velocity distributions discussed here 
with Galactic field stars, rather than the Galactic cluster 
stars discussed in Paper III. Indeed, it has been shown in 
Paper III and by Wolff et al. p007p that the mean rota- 
tional velocities of Galactic cluster objects are higher than 
Galactic field stars. For example, strong peaks in the 
Galactic cluster rotational velocity distributions oc- 
cur at ~200kms~^ while the field population dis- 
cussed below peaks at ^50kms^^. 

Abt et al. PUD^ and Howarth et al. ([TM7|) present ro- 
tational velocities for a large number of B-type and 0-type 
stars in the Galactic field respectively. In order to maintain 
consistency between the Galactic and Magellanic Cloud 
samples we have randomly chosen Galactic objects from 
these catalogues in such a way as to populate an equivalent 
spectral type distribution as the Magellanic Cloud samples, 
with 198 objects being selected in total. In Fig. [10] cumula- 
tive probability functions are plotted for this Galactic sam- 
ple and the LMC and SMC objects shown in Fig. [51 It is 
clear that the SMC objects typically rotate faster than the 
Galactic objects with the LMC distribution lying between 
the two. Given that we are randomly selecting Galactic 
comparison objects, we have carried out the comparison 
1000 times. The Student's t-test and the K-S test give the 
probability that the Galactic and LMC samples are drawn 
from the same distribution to be 43% and 24% respectively 
from the mean of these 1000 comparisons. Similarly the 



two tests show that the probability the Galactic and SMC 
samples are drawn from the same distribution of rotational 
velocities to be 0.3% and 0.8%. Hence while there is some 
evidence that the LMC sample of stars may rotate faster 
that those in the Galactic field, the faster rotation of the 
SMC B-type objects compared to the Galactic objects is 
significant at the 3(t level. 
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Fig. 10. Cumulative probability functions of the rotational 
velocities of the single core-hydrogen burning LMC 
and SMC objects with masses <25Mq. For compari- 
son Galactic (GAL) objects from Howarth et al. (|1997p and 
Abt et al. (|2002p are included (see text). The median in- 
trisic rotational velocities are 125, 135 and 183kms~^ for 
the Galactic, LMC and SMC samples respectively (assum- 
ing random angles of inclination) indicating a metallic- 
ity dependence. 



5.3. Implications for gamma-ray burst progenitors 

The progenitors of long duration gamma-ray bursts (here- 
after GRBs) have been suggested to be rapidly rotating 
massive carbon-oxygen stars that have lost their entire hy- 
drogen and helium envelopes. This poses a problem for stel- 
lar evolutionary models in that one has to remove the outer 
envelope of a massive star, to create a Wolf-Rayet type ob- 
ject (WR), but maintain or enhance the rotational velocity 
(Woosley & Reger ^UU^. 

Very rapidly rotating main-sequence stars have been 
proposed as candidates for chemically homogeneous evo- 
lution, in which the stars are homogeneously mixed on a 
timescale much shorter than their nuclear lifetimes (Maeder 
119871 Langer [r992p . In this case the star becomes a WR star 
directly by evolving to the blue, and its outer hydrogen en- 
velope is effectively mixed into the star rather than being 
lost through stellar winds. This avoids angular momentum 
loss through mass-loss, and hence one can potentially get a 
viable GRB progenitor i.e. a massive CO star that is hydro- 
gen and helium free and rapidly rotating (Yoon & Langer 
i2005t Woosley & Heger 12006)) . Our rotational velocity re- 
sults, and the determination of a frequency distribution can 
provide a critical input parameter for these models. 

Yoon et al. (2006) have recently generated massive star 
evolutionary models at SMC and lower metallicities to 
study the plausibility of chemically homogeneous evolution 
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generating GRB progenitors. They assumed a distribution 
function of rotational velocities which was estimated from 
the most massive 0-type stars in our FLAMES survey data 
set. Mokiem et al. ((200 6| determined projected rotational 
velocities for the 17 unevolved 0-stars in NGC 346 and es- 
timated the intrinsic rotational velocity distributions, as 
in Sect. 15.21 Given the young age of NGC 346 it would 
be expected that angular momentum loss due to winds 
should be small and the current rotational velocity dis- 
tribution should be representative of the initial rotational 
velocity distribution. In order to determine the rate of 
GRB production from rapidly rotating, chemically homo- 
geneous stars, Yoon et al. (.20061^ then assumed that this 
rotational velocity distribution was appropriate to all 
masses and metallicities. They predict that at metallici- 
ties oi Z ^ 0.00001 — 0.001 stars with initial masses as 
low as 11-14 Mq can produce GRBs. By combining our full 
FLAMES survey sample (Mokiem et al., and this paper) we 
have the most extensive rotational velocity survey at low 
metallicities to test these ideas. 



As suggested by Table [9] the more massive stars may ro- 
tate slower than their less massive counterparts, implying a 
possible mass dependence. However, if mass-loss effects are 
small in the 0-star sample whilst on the main-sequence, 
this effect may be due to mass and metallicity dependent 
protostellar winds rather than winds on the main-sequence 
(Yoon et al. 2006). In Fig.[Tl]we plot the cumulative prob- 
ability function for the ratio of projected rotational veloc- 
ity to critical (Keplerian) rotational velocity for the un- 
evolved 0-stars in the Mokiem et al. NGC 346 sample and 
compare this to that derived for the NGC 346 B-type stars 
here. There is a clear offset between the two distributions, 
with the B-stars having higher ratios of rotational to crit- 
ical velocity. We note again that the mass range of the B- 
stars presented here are 10-25 Mq, whereas the 0-stars of 
Mokiem et al. have masses mostly in the range 20-60 Mq. 
This is the first time we have had enough objects to 
carry out this mass-dependent comparison. Hence 
the Yoon et aL rates of GRBs from rapidly rotating 
stars below 20 Mq (greater than wsini/wKcp — 0.5) 
are probably underestimates. 



It would be desirable to recalculate GRB rates 
with our new measurements of lower mass stellar 
rotational velocities. Yoon et al. (in their Fig. 3) cal- 
culate that stars of 30 Mq require a minimum rotational 
velocity of vsini/vkcp ^ 0.27 to be viable GRB candidates, 
which from Fig. [Tl] is about 10-15% of the 0-type stel- 
lar population. A 14 Mq star however requires a minimum 
rotational velocity of vsini/vy^cp 0.43, suggesting about 
15% of that population could be viable GRB progenitors. 
Given the similarity of these percentages, the initial mass 
function should imply more progenitors from the lower mass 
range. We therefore suggest that adopting our ob- 
served rotational velocity for the full mass range 
will increase the fraction of stars which become 
GRBs. A quantatative estimate requires a detailed 
grid calculation as in Yoon et al. This is outside 
the scope of this paper but it should be done with 
the complete rotational velocity distribution that 
we present here. 
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Fig. 11. Cumulative probability functions of projected ro- 
tational velocity to critical rotational velocity ratio for 
the unevolved NGC 346 0-star sample from Mokiem et al. 
P006l Mokiem) and our single core-hydrogen burning 
NGC 346 B-type star sample. This plot shows that the 
rotational velocity distribution of the more massive 
O-stars is peaked towards lower rotational veloci- 
ties than that of the B-type stars. 

6. Conclusions 

Atmospheric parameters and projected rotational velocities 
are presented for approximately 400 O- and early B-type 
stars towards four clusters in the Magellanic Clouds, which 
combines this work, the analysis of some 50 0-type stars 
(Mokiem et al. [MlBl [M?7|) and the analy sis of 100 nar- 
row lined B-type stars (Hunter et al. 120071 and Trundle et 
al. I2007p . This represents the full data set from our ESO 
Large Programme, the VLT Flames survey of massive stars 
in the Magellanic clouds and the Galaxy. From the sub- 
sets of the data we find that supergiants are the 
slowest rotators in the sample, typically having ro- 
tational velocities less than 80kms^. Additionally 
binaries tend to rotate slower than single objects, 
but given that binarity is easier to detect at low 
projected rotational velocities this result may be 
dominated by selection effects. It is also found that 
objects with Be phenomena are typically the fastest 
rotators in our dataset, which is in agreement with 
previous studies. The number of high mass objects 
(M>25Mq) in the combined sample is limited but 
there is some evidence that the most massive ob- 
jects rotate slower than their less massive coun- 
terparts. We utilise the single core-hydrogen burn- 
ing objects along with these subsets to address the 
questions outlined in Sect. [T] 

Current evolutionary models do not predict the 
large number of giants observed in our samples 
which may imply that the theoretical tracks need 
to be extended. Indeed we find that stars classed 
as giants typically rotate some ^50kms^ slower 
than dwarf objects and this is consistent with 
a modest slowing down of the rotational veloc- 
ity over the core hydrogen burning lifetime, in- 
dicating that giants are still in the hydrogen 
burning phase. IVom the bi-model distribution 
of rotational velocities for supergiant and core- 



16 



I. Hunter et al.: Rotational velocity distributions in the LMC and SMC 



Table 13. Peak and width of the Magellanic Cloud rota- 
tional velocity distributions determined by fitting the dis- 
tributions with a Gaussian function, see Sect. \5l2\ 



Region 


Peak 


-!■ half-width 




kms~^ 


kms 


LMC 


100 


150 


SMC 


175 


150 



hydrogen burning objects we observe the end 
of the core-hydrogen burning phase to be at 
3.2 dex. We also show that supergiants are un- 
likely to have evolved directly from the hydrogen 
burning main-sequence, but instead, binarity or 
blue loops may be a more likely explanation for 
the observed supergiant population. 

— Our sample is dominated by field objects and we com- 
pare the rotational velocity distributions with Galactic 
field objects from the literature, finding that the SMC 
metallicity stars rotate fastest, in agreement with the- 
oretical predictions. The difference between the SMC 
stars and Galactic field stars is significant at the 3cr 
level, suggesting we have detected a real difference. 
However we find no significant difference between the ro- 
tational velocity distributions of the Galactic and LMC 
stars. By assuming random angles of inclination, we 
can model the observed v sin i distributions with broad 
Gaussians with the parameters listed in Table fT3l 

— Objects more massive than 25 Mq tend to ro- 
tate slower than less massive objects and this 
could be due to mass-loss effects, with a very massive 
star losing angular momentum through strong stellar 
winds, although this is unexpected in the SMC. An al- 
ternative explanation is that the initial birth rotational 
velocity is mass dependent, perhaps due to protostel- 
lar winds. However, the consequence is that Yoon 
et al. (2006) may have underestimated the ro- 
tational velocity distribution of 10-20 Mq GRB 
progenitors. We therefore suggest that the rotational 
velocity distribution derived here should be employed 
in future GRB rate calculations, and this may even in- 
crease the rate of GRBs expected at low metallicty due 
to the greater prevalence of B-type stars. 
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Table 1. Radial distances (rj-ad, arcmins) of the objects in the NGC346 FLAMES sample with respect to the centre 
of star formation taken as the position of object 435 from the Massey et al. (1989) catalogue of OB stars in NGC346; 
a = 00 59 04.49 6 — -72 10 24.7. The centre of NGC 346 adopted in Paper II for calculation of radial distances was obtained 
from the SIMBAD database (operated at the CDS, Strasbourg, France) and is actually the centre of the ionised shell of 
gas, which is offset from the cluster centre. 
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Table 3. Atmospheric parameters and usini values for the NGC346 sample. In the method column, M06, T06 or H06 
indicates that the values have been directly obtained from Mokiem et al. (|2006[[^07p . Trundle et al. (|2007p or Hunter et 
al. (2007). He indicates the methodology outlined in Sect. 13.21 has been used to derived the atmospheric parameters. 'A' 
indicates that the temperature has been assumed based on spectral type. Values marked with colons indicate parameters 
where fitting has been difficult and the uncertainties may be larger than normal. Where it was possible to measure the 
projected rotational velocity of the secondary object in a double lined spectroscopic binary system, a second entry is 
given with 'B' appended to the object identifier. We do not give atmospheric parameters for those objects listed as Be 
stars (Be), double lined spectroscopic binaries (SB2) or classified with an uncertain spectral type (SP?). 



Star 


Spectral Type 


Tcff 
(K) 


logs 
(dex) 


V sin i 
(kms ^) 


Method 


L/L© 


Mass 
(M©) 


Comments 


NGC 346-001 


07 laf-f 


34100 


3.35 


74 


M06 


6.02 


66 


Binary 


NGC 346-004 


Be(Bl:) 


Be 


- 


266 


- 


5.17 


23 




NGC 346-007 


04 V((f+)) 


42800 


3.95 


120 


M06 


5.45 


40 


Binary 


NGC 346-008 


Ble 


Be 


— 


299 


— 


4.96 


20 




NGC 346-009 


BOe 


Be 


- 


199 


- 


5.08 


24 




NGC 346-010 


07 nin((f)) 


35900 


3.54 


313 


M06 


5.20 


27 
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Bl lb 


24200 


3.20 


<28 


H06 


4.77 


16 
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120 


- 


4.79 


17 


Binary 
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- 


320 


- 






Binary 
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- 


45 


- 


4.79 


17 


Binary 


NGC 346-015B 




SB2 


- 


35 


- 






Binary 
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BO. 5 Vn 
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- 


181 


- 


4.87 


19 


Binary 


NGC 346-018 


09.5 Ille 


32700 


3.33 


138 


M06 


5.10 


24 
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- 


34 


- 


4.69 


16 


Binary 


NGC 346-020B 
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- 


66 


- 






Binary 


NGC 346-021 
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25150 


3.50 


15 


H06 


4.61 


14 




NGC 346-022 


09 V 


36800 


4.20 


55 


M06 


4.95 


24 




NGC 346-023 


B0.2:(Be-Fe) 


Be 


- 


65 


- 


4.81 


18 




NGC 346-024 


B2:shell(Be-Fe) 


Be 


- 


190 


- 


4.59 


19 




NGC 346-025 


09V 


36200 


4.07 


138 


M06 


4.90 


23 


Binary 


NGC 346-026 


BO IV (Nstr) 


32500 


3.75 


75 


He 


4.85 


20 




NGC 346-027 


B0.5 V 


31000 


4.05 


220 


He 


4.78 


18 
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OC6 Vz 


42900 


3.97 


27 


M06 


5.10 


32 




NGC 346-029 


BO V 


32150 


4.10 


25 


H06 


4.82 


19 


Binary 


NGC 346-030 


BO V 


SB2 


- 


183 


- 


4.81 


19 


Binary 


NGC 346-031 


08 Vz 


39500 


3.99 


18 


M06 


4.99 


27 




NGC 346-032 


B0.5 V 


29000 


4.40 


125 


He 


4.68 


17 


Binary 


NGC 346-033 


08 V 


39900 


4.44 


188 


M06 


4.99 


27 




NGC 346-035 


Bl V 


SB2 


- 


145 


- 


4.61 


15 


Binary 


NGC 346-035B 




SB2 


- 


105 


- 






Binary 


NGC 346-036 


B0.5 V(Be-Fe) 


Be 




287 




4.69 


17 




NGC 346-037 


B3 III 


18800 


3.20 


35 


H06 


4.21 


10 




NGC 346-039 


B0.7 V 


25800 


3.60 


20 


H06 


4.51 


13 


Binary 


NGC 346-040 


B0.2 V 


30600 


4.00 


20 


H06 


4.67 


17 


Binary 


NGC 346-041 


B2 (Be-Fe) 


Be 




144 




4.46 


12 




NGC 346-043 


BO V 


33000 


4.25 


10 


H06 


4.71 


18 




NGC 346-044 


Bl II 


23000 


3.50 


40 


H06 


4.33 


11 




NGC 346-045 


B0.5 Vne 


Be 




181 




4.56 


15 




NGC 346-046 


07 Vn 


39700 


4.17 


340 


M06 


4.81 


24 




NGC 346-047 


B2.5 III 


19850 


3.25 


63 


A 


4.15 


10 




NGC 346-048 


Be (B3 shell) 


Be 




158 




4.15 


10 




NGC 346-049 


B8 II 


13000 


2.70 


80 


He 


3.70 


7 




NGC 346-050 


08 Vn 


37200 


4.16 


357 


M06 


4.67 


21 




NGC 346-051 


07 Vz 


41600 


4.33 


18 


M06 


4.87 


27 




NGC 346-052 


B1.5 V 


SB2? 




12 




4.37 


12 


Binary? 


NGC 346-053 


B0.5 V 


29500 


3.75 


170 


He 


4.51 


15 


Binary 


NGC 346-054 


Bl V 


27300 


4.20 


23 


A 


4.41 


13 




NGC 346-055 


B0.5 V 


29500 


4.00 


130 


He 


4.49 


15 




NGC 346-056 


BO V 


31000 


3.80 


15 


H06 


4.55 


16 





20 
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otar 


bpectral iype 


J-eff 

(K) 


logs 
(dex) 


vsini 
(kms ) 


Metnod 


T /T 

L/Lq 


Mass 
(Mo) 


Comments 


NGC 346-057 


B2.5 III 


19850 


3.35 


73 


A 


4.09 


9 




NGC 346-058 


BO. 5 V 


29500 


4.25 


180 


He 


4.47 


14 


Binary 


NGC 346-061 


Bl-2 (Bc-Fe) 


Be 




336 




4.32 


12 




NGC 346-062 


BO. 2 V 


29750 


4.00 


25 


H06 


4.45 


15 




NGC 346-064 


Bl-2 (Bc-Fc) 


Be 




108 




4.31 


12 




NGC 346-065 


B3 (Be-Fe) 


Be 




222 




4.05 


9 




NGC 346-066 


09.5 V 


35600 


4.25 


129 




4.59 


19 




NGC 346-067 


Bl-2 (Be-Fe) 


Be 




351 


~ 


4.29 


12 




NGC 346-068 


BO V (Be-Fe) 


Be 




378 




4.51 


16 




NGC 346-069 


Bl-2 (Be-Fe) 


Be 




186 




4.28 


12 




NGC 346-070 


BO. 5 V 


30500 




109 


[ 1. 


4.43 


15 




NGC 346-072 


Bl-2 (Be-Fe) 


Be 




102 




4.26 


12 




NGC 346-073 


Bl-2 (Be-Fe) 


Be 




190 




4.26 


12 




NGC 346-074 


B3 III 


16500 




52 


1 [■ 


3.80 


7 




NGC 346-075 


Bl V 


27700 


4.30 


10 


H06 


4.31 


13 


Binary 


NGC 346-076 


B2 (Be-Fe) 


Be 




237 




4.21 


11 




NGC 346-077 


09 V 


36500 


3.99 


177 


\lf ifi 


4.65 


20 




NGC 346-078 


B2 III 


SB2? 




154 


A 


4.05 


9 


Binary? 


NGC 346-079 


BO. 5 Vn 


29500 


- ' 


293 


He 


4.37 


14 




NGC 346-080 


Bl V 


27300 


4.25 


216 


A 


4.29 


12 




NGC 346-081 


B2 Illn 


21200 


3.50 


255 


A 


4.03 


9 




NGC 346-082 


B2 III 


21200 


3.70 


168 


A 


4.03 


9 


Binary 


NGC 346-083 


Bl V 


27300 


4.05 


207 


A 


4.28 


12 




NGC 346-084 


Bl V 


27300 


4.25: 


105 


A 


4.27 


12 




NGC 346-085 


B2 III 


SB2 




26 




4.03 


9 


Binary 


NGC 346-085B 




SB2 




11 








Binary 


NGC 346-088 


Bl V 


27300 


: 1 ' 1 : 


84 


\ 


: ;"i 


1 




NGC 346-089 


Bl-2 (Bc-Fc) 


Be 




79 




4.21 


11 




NGC 346-090 


09.5 V 


34900 


4.25 


188 


M06 


4.56 


18 




NGC 346-091 


Blc 


Be 




49 




4.25 


12 


Binary 


NGC 346-092 


Bl Vn 


27300 


3.95 


234 


\ 


4.24 


12 




NGC 346-093 


BO V 


34400 


4.42 


187 


M06 


4.53 


18 




NGC 346-094 


BO. 7 V 


28500 


4.00 


40 


H06 


4.28 


13 




NGC 346-095 


Bl-2 (Be-Fe) 


Be 




227 




4.18 


11 




NGC 346-096 


Bl-2 (Be-Fe) 


Be 




343 




4.18 


11 




NGC 346-097 


09 V 


37500 


4.49 


22 


M06 


4.75 


22 




NGC 346-098 


B1.5 V 


26100 


4.05 


56 


A 


4.16 


11 




NGC 346-099 


B3 III 


18000 


3.40 


94 


He 


3.79 


7 




NGC 346-100 


B1.5 V 


26100 


4.30 


183 


A 


4.16 


11 




NGC 346-101 


Bl V 


27300 


4.25 


29 


A 


4.18 


12 




NGC 346-102 


B3 III 


17700 


3.70 


30 


He 


3.75 


7 




NGC 346-103 


BO. 5 V 


29500 


4.00 


10 


H06 


4.26 


13 




NGC 346-104 


BO V 


33500 


4.45 


17 


He 


4.41 


16 


Binary 


NGC 346-106 


Bl V 


27500 


4.20 


142 


He 


4.17 


12 


Binary 


NGC 346-107 


09.5 V 


35900 


4.23 


55 


M06 


4.40 


18 




NGC 346-108 


B1.5 V 


26100 


3.95 


167 


A 


4.11 


11 




NGC 346-109 


B1.5 V 


26100 


4.25 


123 


A 


4.11 


11 




NGC 346-110 


Bl-2 (Bc-Fe) 


Be 




243 




4.11 


11 


Binary 


NGC 346-112 


09.5 V 


34400 


4.15 


143 


M06 


4.36 


17 




NGC 346-114 


Bl Vn 


27300 


3.90 


287 


A 


4.13 


12 




NGC 346-116 


Bl V 


28250 


4.10 


15 


H06 


4.15 


12 
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Table 4. Atmospheric parameters and usini values of the NGC330 sample. Symbols and labels are the same as given 
in Table O 



Star 


Spectral Type 


Tefl 


log 3 


V sin i 


Method 


L/Lq 


Mass 






(K) 


(dex) 


(kms^^) 






(Mo) 


NGC 330-002 


B3 lb 


14590 


2.15 


<24 


T06 


4.77 


15 


NGC 330-003 


B2 lb 


17210 


2.25 


<36 


T06 


4.89 


16 


NGC 330-004 


B2.5 lb 


17000 


2.30 


<41 


T06 


4.81 


15 


NGC 330-005 


B5 lb 


13700 


2.25 


20 


T06 


4.58 


13 


NGC 330-009 


B5 lb 


13940 


2.45 


28 


T06 


4.45 


12 


NGC 330-010 


B5 lb 


14820 


2.60 


<15 


T06 


4.44 


12 


NGC 330-013 


08.5 III-II((f)) 


34500 


3.40 


78 


M06 


5.40 


32 


NGC 330-014 


B1.5 lb 


20130 


2.75 


45 


T06 


4.69 


14 


NGC 330-016 


B5: II 


14220 


2.60 


40 


T06 


4.24 


10 


NGC 330-017 


B2 II 


22000 


3.35 


21 


T06 


4.66 


14 


NGC 330-018 


B3 II 


18000 


2.95 


48 


T06 


4.45 


12 


NGC 330-020 


B3 II 


16720 


2.85 


<42 


T06 


4.35 


11 


NGC 330-021 


B0.2 III 


30500 


3.70 


204 


He 


4.95 


21 


NGC 330-022 


B3 II 


18860 


3.00 


24 


T06 


4.42 


12 


NGC 330-024 


B5 lb 


14000 


2.70 


46 


He 


4.05 


9 


NGC 330-025 


B1.5e 


Be 


- 


178 


- 


4.64 


15 


NGC 330-026 


B2.5 II 


22500 


3.40 


73 


T06 


4.50 


12 


NGC 330-027 


Bl V 


22040 


3.20 


80 


T06 


4.46 


12 


NGC 330-028 


Bl V 


27300 


3.75 


76 


A 


4.63 


15 


NGC 330-029 


B0.2 V(Be-Fe) 


Be 


- 


209 


- 


4.75 


18 


NGC 330-031 


B0.5 V(Be-Fe) 


Be 


- 


178 


- 


4.69 


17 


NGC 330-032 


B0.5 V 


29700 


4.15 


17 


T06 


4.67 


17 


NGC 330-033 


B1.5 V 


26100 


3.90 


105 


A 


4.53 


14 


NGC 330-034 


Bl-2e 


Be 


- 


231 


- 


4.53 


14 


NGC 330-035 


B3 II 


19000 


3.20 


40 


He 


4.22 


10 


NGC 330-036 


B2 II 


21200 


3.25 


39 


A 


4.33 


11 


NGC 330-038 


Bl V 


27300 


3.75 


150 


A 


4.53 


14 


NGC 330-039 


BO V 


33000 


4.25 


120 


He 


4.72 


18 


NGC 330-040 


B2 III 


21200 


3.25 


106 


A 


4.25 


11 


NGC 330-041 


BO V 


32000 


4.15: 


127 


He 


4.64 


17 


NGC 330-042 


B2 II 


25450 


3.75 


27 


T06 


4.38 


12 


NGC 330-043 


BO V 


33000 


4.10: 


250 


He 


4.64 


18 


NGC 330-044 


Bl-2 (Be-Fe) 


Be 


- 


184 


- 


4.37 


12 


NGC 330-045 


B3 III 


18450 


3.25 


133 


A 


4.01 


9 


NGC 330-046 


09.5 V 


34000 


4.25: 


165 


He 


4.64 


18 


NGC 330-047 


Bl V 


26700 


4.00 


28 


T06 


4.35 


12 


NGC 330-048 


B0.5 V 


29000 


4.00 


73 


He 


4.42 


14 


NGC 330-049 


09 V 


35000 


4.15 


40 


He 


4.64 


19 


NGC 330-050 


B3e 


Be 


- 


214 


- 


4.05 


9 


NGC 330-051 


B1.5 V 


26100 


4.05 


273 


A 


4.30 


12 


NGC 330-052 


08.5 Vn 


35700 


3.91 


291 


M06 


4.60 


19 


NGC 330-053 


B0.5 V 


29650 


4.25 


82 


A 


4.42 


14 


NGC 330-054 


B2 (Be-Fe) 


Be 


- 


147 


- 


4.25 


8 


NGC 330-055 


B0.5 V 


31000 


4.10 


219 


He 


4.46 


15 


NGC 330-056 


B2 III 


21200 


3.50 


108 


A 


4.06 


9 


NGC 330-057 


B0.5 V 


29000 


4.15 


104 


He 


4.34 


13 


NGC 330-058 


B3: 


SP? 




263 




3.97 


10 


NGC 330-059 


B3 III 


18450 


3.25 


123 


A 


3.89 


8 


NGC 330-060 


B2.5 (Be-Fe) 


Be 




88 




4.08 


9 


NGC 330-062 


B3e 


Be 




241 




3.95 


8 


NGC 330-063 


Bl-3 


SP? 




199 




4.17 


11 


NGC 330-064 


B3:e 


Be 




269 




3.93 


8 


NGC 330-065 


Bl-3 (Be-Fe) 


Be 




284 




4.13 


11 


NGC 330-066 


B3 III 


18500 


3.40 


126 


He 


3.83 


8 


NGC 330-067 


B2.5 III 


19850 


3.40 


65 


A 


3.90 


8 



Comments 



Binary 



Binary 



Binary 



Binary 
Binary? 
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Table 4. -continued. 



I. Hunter et al.: Rotational velocity distributions in the LMC and SMC 



star 


opectral iype 


J-eff 


logg 


usini 


Metnod 


T /T 

L/Lq 


Mass 






(K) 


(dex) 


(kms ) 




(Mo) 


NGC 330-068 


B1.5 (Be-Fe) 


Be 


— 


23 


— 


4.16 


11 


NGC 330-069 


B3 III 


19000 


3.40 


193 


He 


3.85 


8 


NGC 330-070 


B0.5e 


Be 




348 


— 


4.29 


13 


NGC 330-071 


B3 III 


19000 




92 


He 


3.84 


8 


NGC 330-072 


BO. 5 V 


29650 


- 


99 


A 


4.29 


13 


NGC 330-073 


B8 lb 


12000 


2.20 


49 


He 


3.37 


5 


NGC 330-074 


BO V 


32020 


4.20 


29 


T06 


4.35 


15 


NGC 330-075 


B8 II 


13500 


3.20 


166: 


He 


3.46 


6 


NGC 330-076 


B3 (Be-Fe) 


Be 




62 


~ 


3.86 


8 


NGC 330-079 


B3 III 


19500 


:;. If 


146 


He 


3.82 


8 


NGC 330-080 


Bl-3 


SP? 




283 




4.06 


10 


NGC 330-081 


Bl-3 


19000 


3.50 


190 


He 


3.79 


7 


NGC 330-082 


Bl-3 


SP? 




192 




4.05 


10 


NGC 330-083 


B3 III 


18000 


3.25 


140 


He 


3.72 


7 


NGC 330-084 


B3 V-III 


SB2? 




26 




3.83 


8 


NGC 330-085 


B3:e 


Be 




191 




3.82 


8 


NGC 330-086 


B2.5 III 


19850 


3.50 


89 


A 


3.80 


8 


NGC 330-087 


Be-Fe 


Be 




214 


~ 


3.81 


8 


NGC 330-089 


Bl-5 


18500 




193 


He 


3.70 


7 


NGC 330-090 


B3 III 


18450 


3.30 


167 


A 


3.69 


7 


NGC 330-091 


BOc 


Be 




272: 




4.25 


14 


NGC 330-094 


Bl-5 


SP? 




330 




3.77 


8 


NGC 330-095 


B3 III 


18450 


3.45 


20 


A 


3.67 


7 


NGC 330-096 


Bl-3 (Be-Fe) 


Be 




175: 




3.97 


10 


NGC 330-097 


Bl V 


27300 


.).!)' 


154 


A 


4.05 


11 


NGC 330-098 


B0.2: V 


SP? 




148 




4.18 


13 


NGC 330-099 


B2-3 


SP? 




187 




3.86 


9 


NGC 330-100 


Be (BO-3) 


Be 




373 




3.99 


10 


NGC 330-101 


B2.5 III 


19850 


3.50 


48 


A 


3.73 


7 


NGC 330-102 


B2-3 III 


19850 


3.60 


129 


A 


3.72 


7 


NGC 330-103 


Bl-3 


SP? 




131 




3.94 


9 


NGC 330-104 


BO: V 


SP? 




309: 




4.19 


14 


NGC 330-105 


Bl-3 


SP? 




13 




3.91 


9 


NGC 330-106 


Bl-2 


SP? 




71 




3.95 


10 


NGC 330-107 


B3: V-III 


16500 


3.25 


34 


He 


3.49 


6 


NGC 330-108 


B5 III 


14500 


3.30 


102 


He 


3.35 


5 


NGC 330-109 


B3 III 


18450 


3.20 


213 


A 


3.60 


7 


NGC 330-110 


B2 III 


21200 


3.40 


110 


A 


3.74 


8 


NGC 330-111 


Bl-3 


SP? 




173 




3.79 


8 


NGC 330-112 


Bl-3 (Be-Fe) 


Be 




262 




3.79 


8 


NGC 330-113 


Bl-3 


SP? 


— 


327 


— 


3.78 


8 


NGC 330-114 


B2 III 


23800 


3.90 


17 


T06 


3.83 


9 


NGC 330-116 


B3 III 


19500 


3.55 


4 


He 


3.62 


7 


NGC 330-118 


Bl-2 


SP? 




32 




3.89 


10 


NGC 330-119 


Bl-3 


SP? 




265 




3.84 


9 


NGC 330-120 


B3: V-III 


18500 


3.40 


137: 


He 


3.47 


6 


NGC 330-123 


09.5 V 


35000 


4.25: 


26 


He 


4.58 


18 


NGC 330-124 


B0.2 V 


30980 


4.25 


95 


T06 


4.42 


15 


NGC 330-125 


B2 III 


21200 


3.50 


86 


A 


4.01 


9 



Comments 



Binary 



Binary 



Binary 
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Table 5. Atmospheric parameters and wsini values of the Nil sample. Symbols and labels are the same as given in 
Table El 



Star 


Spectral Type 


Tcff 




V sin i 


Method 


— T — 7t 

L/Lq 


Mass 


Comments 






(K) 


(dex) 


(kms^^) 




(Mo) 




N 11-001 


B2 la 


18750 


2.50 


<52 


H06 


5.66 


38 




N 11-002 


B3 la 


15800 


2.10 


<53 


H06 


5.26 


24 




N 11-003 


Bl la 


23200 


2.75 


<54 


H06 


5.42 


30 




N 11-004 


09.7 lb 


31600 


3.36 


40 


M06 


5.80 


48 




N 11-008 


BO. 5 la 


25450 


3.00 


<43 


H06 


5.39 


30 




N 11-009 


B3 lab 


15000 


2.15 


30 


H06 


4.85 


17 


Binary 


N 11-010 


09.5 III -f Bl-2 


SB2 


- 


152 


- 


5.59 


40 


Binary 


Nll-OlOB 




SB2 


- 


12 


- 


- 


- 


Binary 


N 11-011 


OC9.5 II 


29500 


3.15 


<60 


He 


5.49 


35 


Binary 


N 11-012 


Bl la 


20500 


2.55 


40 


H06 


5.13 


22 




N 11-014 


B2 lab 


19100 


2.55 


30 


H06 


5.03 


19 


Binary 


N 11-015 


BO. 7 lb 


23600 


2.95 


58 


H06 


5.23 


24 




N 11-016 


Bl lb 


21700 


2.75 


54 


H06 


5.13 


22 




N 11-017 


B2.5 lab 


16500 


2.30 


<25 


H06 


4.82 


17 


Binary 


N 11-023 


BO. 7 lb 


24000 


2.90 


66 


H06 


5.09 


21 




N 11-024 


Bl lb 


21600 


2.80 


<30 


H06 


4.96 


18 




N 11-026 


02.5 III(f*) 


53300 


4.00 


109 


M06 


5.92 


82 




N 11-029 


OC9.7 lb 


28750 


3.30 


40 


H06 


5.21 


25 




N 11-031 


0N2 III(f*) 


45000 


3.85 


116 


M06 


5.84 


61 




N 11-032 


07 11(f) 


35200 


3.45 


<65 


M06 


5.43 


34 




N 11-033 


BO Illn 


27200 


3.21 


256 


M06 


5.07 


21 




N 11-034 


BO. 5 III 


25500 


3.25 


203 


He 


5.03 


20 


Binary 


N 11-035 


09 11(f) 


31000 


3.25 


<56 


He 


5.23 


27 


Binary 


N 11-036 


BO. 5 lb 


23750 


3.10 


30 


H06 


4.95 


18 




N 11-037 


BO III 


28100 


3.25 


100 


H06 


5.08 


23 


Binary 


N 11-038 


05 III(f+) 


41000 


3.72 


145 


M06 


5.69 


48 




N 11-039 


B2 III 


21700 


3.00 


157 


A 


4.81 


16 


Binary? 


N 11-040 


BO: Illn 


29500 


3.30 


273 


He 


5.11 


24 




N 11-042 


BO III 


29000 


3.60 


30 


H06 


5.05 


22 


Binary 


N 11-045 


09-9.5 III 


32300 


3.32 


105 


M06 


5.15 


25 




N 11-046 


09.5 V 


33500 


4.25 


205 


He 


5.21 


28 


Binary 


N 11-047 


BO III 


29200 


3.65 


55 


H06 


5.03 


22 


Binary 


N 11-048 


06.5 V((f)) 


40700 


4.19 


130 


M06 


5.38 


37 




N 11-051 


05 Vn((f)) 


42400 


3.75 


333 


M06 


5.31 


36 




N 11-052 


09.5 V 


SB2 




16 


- 


5.23 


29 


Binary 


N 11-054 


Bl lb 


23500 


3.05 


45 


H06 


4.79 


16 




N 11-056 


Ble 


Be/SB2 


- 


205 


- 


4.89 


19 


Binary 


N11-056B 




Be/SB2 


- 


25 


- 


- 


- 


Binary 


N 11-058 


05.5 V((f)) 


41300 


3.89 


85 


M06 


5.27 


34 




N 11-059 


09 V 


34000 


3.85 


103 


He 


5.13 


26 


Binary 


N 11-060 


03 V ((f*)) 


45700 


3.92 


106 


M06 


5.57 


49 




N 11-061 


09 V 


33600 


3.51 


87 


M06 


5.20 


27 




N 11-062 


B0.2 V 


30400 


4.05 


25 


H06 


4.95 


21 


Binary 


N 11-063 


09: V 


35000 


4.25 


196 


He 


5.11 


26 




N 11-064 


BO. 2: Vn 


SB2? 


- 


206 


- 


4.92 


20 


Binary? 


N 11-065 


06.5 V((f)) 


41700 


3.89 


83 


M06 


5.17 


33 




N 11-066 


07 V((f)) 


39300 


3.87 


71 


M06 


5.10 


30 




N 11-068 


07 V((f)) 


39900 


4.13 


54 


M06 


5.06 


29 




N 11-069 


Bl III 


24300 


3.30 


80 


H06 


4.63 


15 




N 11-070 


B3 III 


19500 


3.30 


62 


He 


4.40 


12 




N 11-072 


B0.2 V 


28800 


3.75 


15 


H06 


4.77 


18 




N 11-073 


BO. 5 (Be-Fe) 


Be 




222 




4.78 


18 




N 11-074 


BO. 5 (Be-Fe) 


Be 




169 




4.78 


18 




N 11-075 


B2 III 


21800 


3.35 


25 


H06 


4.48 


12 


Binary 


N 11-076 


BO. 2 la 


26500 


2.90 


72 


He 


4.67 


16 
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Table 5. -continued. 
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otar 


Spectral Type 


J-eff 


logs 


V sin i 


Method 


T /T 


Mass 


Comments 






(K) 


(dex) 


(kms ) 




(Mo) 




N 11-077 


B2 III 


21500 


3.30 


117 


He 


4.47 


12 


Binary 


N 11-078 


B2 (Bc-Fc) 


Be 




93 




4.59 


14 


Binary 


N 11-079 


BO. 2 V 


32500 


4.25 


38 


[ 1. 


4.88 


21 




N 11-081 


B0:n (Be-Fe) 


Be 




363 


He 


4.79 


19 




N 11-082 


Bl-2 -l-early B 


SB2 




95 




4.58 


14 


Binary 


N11-082B 




SB2 




110 




— 




Binary 


N 11-083 


B0.5 V 


29300 




20 


1 Ii if) 


4.71 


1 7 


Binary 


N 11-084 


B0.5 V 


29500 


4.25 


108 


He 


4.75 


18 


Binary 


N 11-085 


B0.5 V 


SB2 


— 


105 


— 


4.73 


17 


Binary 


N 11-085B 




SB2 




28 




~ 


~ 


Binary 


N 11-086 


Bl V 


26800 


4.25 


75 


A 


4.65 


16 




N 11-087 


09.5 Vn 


32700 


4.04 


276 


M06 


4.91 


21 




N 11-088 


Bl III 


24150 


3.45 


240 


A 


4.54 


14 




N 11-089 


B2 III 


21700 


3.20 


117 


A 


4.42 


12 


Binary 


N 11-090 


B2e 


Be 




104 




4.53 


14 




N 11-093 


B2.5 III 


19500 


3.30 


73 


[ li 


4.29 


11 




N 11-094 


Bl III 


SB2 


— 


147 




4.50 


13 


Binary 


N 11-095 


Bl Vn 


26800 


3.85 


267 


\ 


4.59 


15 




N 11-096 


Bl III 


SB2 




32 




4.49 


13 


Binary 


N 11-097 


B3 II 


17500 


2.95 


52 


[ [■ 


4.17 


10 




N 11-098 


B2 III 


21700 


3.50 


45 


A 


4.37 


12 


Binary 


N 11-100 


B0.5 V 


29700 


4.15 


30 


H06 


4.68 


17 




N 11-101 


BO. 2 V 


29800 


3.95 


70 


H06 


4.68 


17 




N 11-102 


BO. 2 V 


31000 


4.20 


218 


He 


4.73 


18 




N 11-103 


Bl-2 -1- early B 


SB2 




125 




4.53 


14 


Binary 


N 11-103B 




SB2 




80 




~ 




Binary 


N 11-104 


B1.5 V 


25700 


3.75 


153 


\ 


4.52 


1 I 




N 11-105 


Bl V 


SB2? 


— 


116 


— 


4.56 


15 


Binary? 


N 11-106 


BO V 


31200 


4.00 


25 


H06 


4.72 


18 




N 11-107 


Bl-2 -h early B 


SB2 




80 




4.51 


13 


Binary 


N 11-107B 




SB2 




170 




~ 




Binary 


N 11-108 


09.5 V 


32150 


: 


25 


■ ' 


4.73 


' 




N 11-109 


B0.5 lb 


25750 


3.20 


40 


H06 


4.48 


12 




N 11-110 


Bl III 


23100 


3.25 


25 


H06 


4.37 


12 




N 11-111 


B1.5 III 


22950 


3.50 


101 


A 


4.36 


12 




N 11-113 


B0.5 III 


SB2 




38 




4.49 


14 


Binary 


N 11-113B 




SB2 




20 




— 


- 


Binary 


N 11-114 


BO Vn 


32500 


- 


299 


II 


4.70 


19 




N 11-115 


Bl III 


24150 


3.65 


53 


A 


4.39 


12 




N 11-116 


B2 III 


21700 


3.30 


160 


A 


4.28 


11 




"NT 11 1 1 
IN 11-11/ 


DL Vn 


ZDoUU 


/I on 
4.20 


O /I T 


A 

A 


4.4D 


1 o 

13 




N 11-118 


B1.5 V 


25700 


3.80 


150 


A 


4.42 


13 


Binary 


N 11-119 


B1.5 V 


SB2 




259 




4.42 


13 


Binary 


N 11-120 


B0.2 Vn 


31500 


4.10 


207 


He 


4.63 


17 




N 11-121 


Bl Vn 


26800 


4.20 


265 


A 


4.45 


14 




N 11-122 


09.5 V 


33000 


4.10 


173 


He 


4.62 


18 




N 11-123 


09.5 V 


34800 


4.22 


110 


M06 


4.58 


19 




N 11-124 


B0.5 V 


28500 


4.20 


45 


H06 


4.47 


14 


Binary 
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Table 6. Atmospheric parameters and vs'mi values of the NGC 2004 sample. Symbols and labels are the same as given 
in Table O 



Star Spectral Type T^f logg v sin i Method L/Lq Mass Comments 







(K) 


(dex) 


(kms^i) 






(Mo) 




NGC 2004-003 


B5 la 


14450 


2.10 


33 


T06 


5.04 


20 


Binary 


NGC 2004-005 


B8 la 


12390 


1.90 


30 


T06 


4.95 


18 


Binary 


NGC 2004-007 


B8 Ila 


12250 


2.00 


<30 


T06 


4.90 


17 




NGC 2004-010 


B2.5 lab 


17160 


2.40 


<30 


T06 


5.02 


19 




NGC 2004-011 


B1.5 la 


21250 


2.75 


<40 


T06 


5.22 


24 




NGC 2004-012 


B1.5 lab 


21270 


2.87 


<40 


T06 


4.92 


18 




NGC 2004-013 


B2 II 


21700 


3.20 


145 


A 


4.94 


18 




NGC 2004-014 


B3 lb 


17800 


2.85 


<30 


T06 


4.72 


15 




NGC 2004-015 


B1.5 II 


22950 


3.15 


165 


A 


4.96 


18 


Binary 


NGC 2004-020 


B1.5 II 


22950 


3.20 


145 


A 


4.91 


18 


Binary 


NGC 2004-021 


B1.5 lb 


21450 


3.00 


<40 


T06 


4.82 


16 




NGC 2004-022 


B1.5 lb 


21780 


3.15 


<30 


T06 


4.79 


16 




NGC 2004-023 


B2 (Be-Fe) 


Be 


- 


102: 


- 


4.85 


17 


Binary? 


NGC 2004-024 


B1.5 Illn 


22950 


3.05 


240 


A 


4.72 


15 




NGC 2004-025 


B2 (Be-Fe) 


Be 


- 


83 


- 


4.75 


15 


Binary 


NGC 2004-026 


B2 II 


22900 


3.65 


<20 


T06 


4.68 


15 


Binary 


NGC 2004-027 


BOe 


Be 


- 


182: 


- 


5.00 


23 




NGC 2004-029 


B1.5e 


23100 


3.50 


30 


T06 


4.65 


14 


Binary 


NGC 2004-030 


B0.2 lb 


29000 


3.75 


123 


He 


4.87 


19 


Binary 


NGC 2004-031 


B2 II 


21700 


3.35 


65 


A 


4.58 


13 


Binary 


NGC 2004-032 


B2 II 


21700 


3.40 


110 


A 


4.58 


13 




NGC 2004-035 


Bl: (Be-Fe) 


Be 


- 


244: 


- 


4.74 


17 




NGC 2004-036 


B1.5 III 


22870 


3.35 


42 


T06 


4.58 


13 




NGC 2004-039 


B1.5e 


Be 


- 


212: 


- 


4.67 


15 




NGC 2004-041 


B2.5 III 


20450 


3.30 


101 


A 


4.43 


12 


Binary 


NGC 2004-042 


B2.5 III 


20980 


3.45 


42 


T06 


4.45 


12 




NGC 2004-043 


B1.5 III 


22950 


3.50 


24 


A 


4.52 


13 




NGC 2004-044 


B1.5: 


SB2 


- 


18 


- 


4.62 


15 


Binary 


NGC 2004-044B 




SB2 


- 


55 


- 


- 


- 


Binary 


NGC 2004-045 


B2 III 


21700 


3.50 


128 


A 


4.46 


12 


Binary 


NGC 2004-046 


B1.5 III 


26090 


3.85 


32 


T06 


4.62 


15 




NGC 2004-047 


B2 III 


21700 


3.35: 


133 


A 


4.42 


12 


Binary 


NGC 2004-048 


B2.5e 


Be 


- 


173 


- 


4.49 


13 




NGC 2004-050 


B2.5 III 


20450 


3.30 


109 


A 


4.35 


11 


Binary 


NGC 2004-051 


B2 III 


21700 


3.40 


70 


A 


4.40 


12 




NGC 2004-052 


B2 III 


21700 


3.60 


138 


A 


4.40 


12 




NGC 2004-053 


B0.2 Ve 


31500 


4.15 


7 


T06 


4.77 


18 




NGC 2004-054 


B2 III 


21700 


3.40 


114 


A 


4.39 


12 


Binary 


NGC 2004-055 


B2.5 III 


20450 


3.30 


118 


A 


4.33 


11 




NGC 2004-056 


B1.5e 


Be 


- 


229 


- 


4.52 


14 




NGC 2004-058 


09.5 V (Nstr) 


33500 


4.10 


195 


He 


4.80 


20 


Binary 


NGC 2004-059 


B2 III 


21700 


3.45 


91 


A 


4.35 


11 


Binary 


NGC 2004-060 


B2 III 


21700 


3.45 


134 


A 


4.35 


11 




NGC 2004-061 


B2 III 


20990 


3.35 


40 


T06 


4.31 


11 




NGC 2004-062 


B0.2 V 


30250 


4.35 


106 


A 


4.67 


17 




NGC 2004-063 


B2 III 


21700 


3.50 


107 


A 


4.32 


11 




NGC 2004-064 


B0.7-B1 III 


25900 


3.70 


28 


T06 


4.48 


13 




NGC 2004-065 


B2.5 III 


20450 


3.40 


223 


A 


4.25 


11 


Binary? 


NGC 2004-066 


B1.5 Vn 


25700 


3.70 


238 


A 


4.48 


13 




NGC 2004-067 


B1.5e 


Be 




237 




4.48 


13 




NGC 2004-068 


B2.5 III 


20450 


3.65 


62 


A 


4.25 


11 




NGC 2004-069 


B0.7 V 


28000 


4.00 


178 


He 


4.54 


15 




NGC 2004-070 


B0.7-B1 III 


27400 


3.90 


46 


T06 


4.51 


14 




NGC 2004-071 


B1.5 III 


22950 


3.50 


98 


A 


4.33 


12 




NGC 2004-073 


B2 III 


21700 


3.55 


37 


A 


4.26 


11 
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Table 6. -continued. 



I. Hunter et al.: Rotational velocity distributions in the LMC and SMC 



btar 


bpectral iype 


J-eff 


logs 


vsini 


Metnoa 


T If 

L/Lq 


Mass 


Comment 






(K) 


(dex) 


(kms ) 




(Mo) 




NGC 2004-074 


B0.7-B1 V 


27375 


4.25 


130 


A 


4.49 


14 


Binary 


NGC 2004-075 


B2 III 


21700 


3.55 


116 


A 


4.26 


11 




NGC 2004-076 


B2.5 III 


20450 


3.65 


37 


A 


4.20 


10 




NGC 2004-077 


BO. 5 V 


29500 


4.00: 


215 


He 


4.56 


16 




NGC 2004-078 


B2 III 


21700 


3.55 


115 


A 


4.26 


11 


Binary 


NGC 2004-079 


B2 III 


21700 


3.60 


165 


A 


4.25 


11 


Binary 


NGC 2004-080 


B2.5 III 


20450 


3.40 


85 


A 


4.19 


10 




NGC 2004-081 


Bl V 


26800 


4.00 


105: 


A 


4.45 


13 




NGC 2004-082 


B1.5 V 


25700 


4.10 


161 


A 


4.42 


13 




NGC 2004-083 


B1.5:c 


Be 




194 




4.42 


13 


Binary 


NGC 2004-084 


B1.5 III 


27395 


4.00 


36 


1 : ' 


4.46 


14 




NGC 2004-085 


B2.5 III 


20450 


3.40 


150 


A 


4.18 


10 




NGC 2004-086 


B2 III 


21700 


3.85 


14 


A 


4.24 


11 




NGC 2004-087 


B1.5 V 


25700 


4.40 


35 


A 


4.40 


13 




NGC 2004-088 


B2.5 III 


20450 


3.65 


200: 


A 


4.18 


10 


Binary 


NGC 2004-089 


B2.5e 


Be 




288 




4.31 


12 




NGC 2004-090 


09.5 III 


32500 


1.1) 


16 


1 lifi 


4.64 


17 




NGC 2004-091 


B1.5 III 


26520 


4.05 


40 


T06 


4.42 


13 


Binary 


NGC 2004-092 


B2e 


Be 


— 


171 


— 


4.34 


12 




NGC 2004-093 


B3 III 


20000 


3.65: 


189: 


He 


4.14 


10 


Binary 


NGC 2004-094 


B2.5 III 


20450 


3.40 


84: 


A 


4.16 


10 


Binary 


NGC 2004-095 


B1.5 V 


25700 


4.10 


138 


A 


4.38 


13 




NGC 2004-096 


B1.5e 


Be 




245 




4.38 


13 




NGC 2004-097 


B2 III 


21700 


3.60 


62 


\ 


4.20 


10 




NGC 2004-098 


B2 III 


21700 


3.75 


90 


A 


4.20 


10 




NGC 2004-099 


B2 III 


21700 


3.40 


119 


A 


4.19 


10 




NGC 2004-100 


Bl Vn 


26800 


3.70 


323 


A 


4.38 


13 




NGC 2004-101 


B2 III 


21700 


3.45 


131 


A 


4.18 


10 




NGC 2004-102 


B2 III 


SB2? 


~ 


64: 




4.18 


10 


Binary 


NGC 2004-103 


B2 III 


21700 


3.85 


35 


\ 


4.18 


10 




NGC 2004-104 


B1.5 V 


25700 


3.90 


274 


A 


4.34 


12 




NGC 2004-105 


B1.5 V 


25700 


3.90 


235 


A 


4.34 


12 




NGC 2004-106 


B2 III 


21700 


3.50 


41 


A 


4.17 


10 




NGC 2004-107 


BO. 5 V 


28500 


3.90 


146 


He 


4.43 


14 


Binary 


NGC 2004-108 


B2.5 III 


22600 


4.00 


13 


T06 


4.21 


10 


Binary 


NGC 2004-109 


B2.5 III 


20450 


3.50 


41 


A 


4.11 


10 


Binary 


NGC 2004-110 


B2 III 


21700 


3.40 


121 


A 


4.17 


10 




NGC 2004-111 


B2.5 III 


20450 


3.30 


55 


A 


4.08 


9 




NGC 2004-112 


B2 III 


21700 


3.65 


72 


A 


4.14 


10 




NGC 2004-113 


B2.5 Illn 


20450 


3.25 


307: 


A 


4.08 


9 




NGC 2004-114 


B2 III 


21700 


3.60 


59 


A 


4.12 


10 




NGC 2004-115 


B2e 


Be 




8 




4.24 


11 


Binary 


NGC 2004-116 


B2 III 


21700 


3.55 


63 


A 


4.12 


10 




NGC 2004-117 


B2 III 


21700 


3.60 


75 


A 


4.11 


10 




NGC 2004-118 


B1.5 V 


SB2? 




90: 




4.25 


12 


Binary 


NGC 2004-119 


B2 III 


23210 


3.75 


15 


T06 


4.15 


10 





